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Chapter  1 

INTRODUCTION 


Vibration  control  has  always  been  an  issue  of  concern  for  engineers.  When 
one  thinks  of  anything  which  is  engineered,  one  generally  encounters  vibrations. 
Vibrations  can  be  a  problem  for  buildings  that  house  rotating  machinery,  buildings 
that  withstand  wind  loads,  computers,  and  road  surfaces.  This  study  deals, 
specifically,  with  active  control  of  structural  vibration  intensity  in  a  beam  and  strives 
to  deal  with  the  larger  problem  of  active  vibration  control. 

Section  1.1  gives  a  brief  history  of  adaptive  control  of  vibrations.  Section  1.2 
gives  a  history  of  measurement  and  adaptive  control  of  vibration  intensity.  Section 
1.3  explains  the  worir  that  was  performed  and  the  goals  of  the  study.  Section  1.4 
details  the  structure  of  the  thesis. 

1.1  Hbtonr  of  Adaptive  Control 

Adaptive  control  can  be  traced  t»ck  as  for  as  the  1950’s  when  "Drenick  and 
Shahbender  borrowed  the  term  'adaptive'  foom  biology  where  it  describes  the  ability 
of  an  organism  to  adjust  itself  to  its  environment”  [1,2].  In  1958,  Anderson  et  al 
[3]  decided  that  an  adaptive  system  must  perform  "a  continuous  measurement  of 
system  dynamic  performance...,  a  continuous  evaluation  of  the  dynamic  performance 
on  the  basis  of  some  i»edetermined  criterion,"  and  "a  continuous  readjustment  of 
system  control  parameters  for  'optimum'  operation...'’  Furthermore,  in  1958,  Grogin- 
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sky  [4]  decided  that  an  adaptive  system  should  outperform  a  conventional  system  by 
"allowing  modifications  of  the  control  system  while  the  system  is  in  operation." 

Non-adaptive  active  control  of  structural  vibrations  was  studied  in  the  1950’s, 
when  Olson  [5]  generalized  his  "electronic  sound  absorber"  to  control  vibrations. 
He  called  his  device  an  "electronic  vibration  reducer,"  and  he  stated  that  it  could  "be 
used  to  isolate  vibrating  machines..."  Adaptive  control  of  structural  vibrations  was 
studied  in  the  1960’s;  Booesho  and  Bollinger  [6]  described  how  to  build  a  "self- 
optimizing  vibration  damper."  The  device  used  an  analog  control  circuit  to  vary  the 
damper’s  stiffness.  The  control  criterion  desired  to  have  the  "phase  angle  6  between 
motion  of  the  damper  mass  m  and  main-system  mass  M  ...  [be]  90°..." 

With  recent  developments  in  computer  technology,  most  notably  speed  of 
processing,  active  control  studies  are  becoming  more  prevalent.  Studies  in  the  1980’s 
investigated  problems  of  vibrations  from  systems  as  simple  as  a  beam  to  those  as 
complex  as  a  "flexible  space  telescope"  [7].  Many  control  theories  were  and  are 
being  examined.  Modal  control  and  distributed-parameter  control  are  replacing  the 
single  parameter  control  schemes  of  the  19S0’s  [8,9,10]. 

In  recent  years  much  research  has  been  performed  generalizing  the  filtered-x 
algorithm  developed  by  Widrow  and  Steams  [11:  288-294],  In  1990,  Sommerfeldt 
and  Tichy  [12]  developed  an  algorithm  based  on  this  filtered-x  algorithm.  Their 
control  approach  provided  for  real-time  system  identification  and  for  multiple  error, 
multiple  output  control.  In  1991,  Sommerfeldt  [13]  applied  this  algorithm,  most 
notably,  to  the  case  of  4  error  sensors  and  4  control  actuators.  A  similar  algorithm 
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has  also  been  developed  to  control  sound  fields;  Sommerfeldt  and  Nashif  [14],  in 
1992,  used  a  filtered-x  algorithm  to  investigate  the  differences  between  controlling  the 
total  acoustic  energy  density  and  controlling  the  squared  pressure. 

1.2  Vibration  Intensity 

This  section  describes  the  history  of  studies  dealing  with  the  measurement  or 
control  of  vibration  intensity.  Section  1.2.1  gives  a  brief  deflnition  of  terms  used  in 
this  thesis.  Section  1.2.2  gives  a  summary  of  previous  studies. 

1.2.1  Deflnition  of  Terms 

Hie  first  term  which  needs  to  be  deHned  is  structural  vibration  intensity. 
Vibration  intensity  in  a  beam  could  also  be  referred  to  as  mechanical  power.  Any 
power  (electrical,  mechanical,  etc.)  per  unit  area  is  referred  to  as  intensity.  Due  to 
the  constant  area  of  a  beam,  having  power  flow  in  one  dimension,  the  power  and 
intensity  vary  by  this  constant  area.  Therefore,  that  which  is  referred  to  as  intensity 
in  this  thesis  is  actually  power. 

The  total,  instantaneous  intensity  at  a  given  point  is  a  sum  of  two  terms:  the 
force  intensity  and  the  moment  intensity.  The  force  intensity  is  Fv-where  F  is  the 
force,  and  v  is  the  velocity-and  the  moment  intensity  is  AfQ-where  M  is  the  moment, 
and  Q  is  the  angular  velocity.  The  total,  instantaneous  intensity  is  represented  as 
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n=Fv+A/Q.  (1-1) 

The  instantaneous  intensity  can  be  divided  between  active  and  reactive 
intensity.  The  active  intensity  is  that  which  is  propagating  down,  and  out  of,  the 
beam  (assuming  there  is  any  propagating  intensity).  The  active  intensity  results  from 
the  force  and  velocity-and  the  moment  and  angular  velocity-being  in  phase,  and  it 
has  a  time-averaged  dc  value.  The  reactive  intensity  is  the  intensity  due  to  any 
standing  waves  in  the  beam.  The  time-averaged  reactive  intensity  has  a  frequency 
twice  that  of  the  excitation  vibration.  Because  the  reactive  intensity  is  periodic  in 
time,  it  could  be  eliminated  by  taking  a  time  average  of  the  total  intensity. 

The  time-averaged  intensity  is  another  quantity  that  will  be  studied.  It  can  be 
expressed  as 

<n\^ljte{FVUljte{MQ^).  (1.2) 

2  2 


1.2.2  Literature  Concerning  Vibration  Intensity 

There  have  been  several  studies  of  measurements  of  vibration  intensity.  In 
1970,  Notseux  [15]  described  a  method  for  measuring  the  "power  flow"  in  beams. 
He  showed  that,  "at  least  for  simple  types  of  waves,"  the  two  terms  of  the  intensity, 
force  and  moment,  are  equal  in  the  far  field,  and  he  showed  that  the  moment  term 
could  be  measured  using  rotational  accelerometers  and  strain  gauges.  In  1984,  Pavic 
[16]  described  a  method  for  measuring  the  intensity  using  an  array  of  velocity 
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transducers.  His  scheme  employed  a  four-point  finite  difference  method  to  calculate 
the  spatial  derivatives  needed  to  calculate  structural  intensity.  Most  recently,  Suen 
[17]  studied  several  finite  difference  methods  that  could  be  used  to  compute  these 
spatial  derivatives.  He  used  a  scanning  laser  vibrometer  to  measure  the  velocity  along 
a  beam,  and  then  he  calculated  time-averaged  intensity. 

Fewer  studies  have  investigated  the  control  of  vibration  intensity.  Redman- 
White  et  al.  [18]  investigated  the  "active  control  of  flexural  wave  power  flow"  in 
1987,  and  concluded  "two  closely  spaced  secondary  force  inputs  can  be  used  to 
absorb  energy  associated  with  the  propagation  of  flexural  waves  in  a  uniform  thin 
beam."  In  1991,  Pan  and  Hansen  [19]  performed  an  analysis  of  the  control  of 
power  flow.  Their  study  included  an  investigation  into  the  effect  of  the  orientation 
of  the  control  actuator.  Numerical  analysis  of  a  beam  structure  was  performed  with 
a  variation  of  the  angle  of  impact  of  the  control  force.  Their  results  showed,  notably, 
that,  in  order  to  control  power  flow  effectively,  the  control  force  must  be  applied 
within  a  few  degrees  of  the  normal. 

1 J  Goals  of  Thesis  Research 

The  goal  of  this  thesis  is  to  develop  an  algorithm  through  which  the  total, 
instantaneous  intensity  can  be  measured  and  controlled.  The  study  investigates 
methods  for  measuring  the  instantaneous  intensity  and  develops  an  appropriate  control 
algorithm.  The  study  also  deals  with  implementing  the  measurement  and  control 
algorithms  and  the  problems  associated  with  this  implementation. 
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The  ultimate  goal  of  this  study  is  to  control  the  active  intensity  in  a  general 
structure.  However,  initially,  it  was  decided  to  control  the  total  intensity,  both  active 
and  reactive.  Also  the  structure  used  is  a  beam  so  as  to  limit  the  control  study  to  one 
dimension. 

1.4  Oreanization  of  the  Thesis 

Chapter  2  develops  a  more  complete  understanding  of  the  intensity.  The 
partial  derivatives  needed  to  calculate  the  intensity  are  explored,  and  the  question  of 
why  one  might  want  to  control  intensity  is  addressed.  Chapter  3  investigates  the 
terms  necessary  for  calculating  the  intensity.  It  describes  the  finite  difference 
schemes  in  both  space  and  time.  Also  considered  in  Chapter  3  is  the  necessity  of 
calibrating  the  accelerometer  array,  as  well  as  additional  calculation  strategies  that 
were  investigated.  Chapter  4  describes  the  development  of  the  control  algorithms. 
The  projection  algorithm,  used  for  system  identirication  is  described  rirst.  The 
filtered-x  algorithm,  used  for  control,  is  described  in  general.  Then  the  specific 
control  schemes  for  implementing  the  riitered-x  algorithm  are  developed.  Chapter  5 
gives  a  description  of  the  experimental  method.  A  description  of  the  experimental 
hardware  is  given  first  This  is  followed  by  a  description  of  what  data  is  collected 
and  in  what  manner  it  is  collected.  Chapter  6  is  a  discussion  of  the  experimental 
results.  The  validity  of  the  intensity  measurement  is  discussed,  and  the  results  of  the 
various  control  strategies  are  also  discussed.  Chapter  7  deals  with  conclusions  drawn 
from  the  study.  Several  recommendations  for  further  research  are  also  discussed. 


Chapter  2 

DESCRIPTION  OF  INTENSITY 


This  chapter  develops  a  more  complete  understanding  of  the  structural 
intensity.  In  Section  2.1,  expressions  for  the  total,  instantaneous  intensity  and  the 
time-averaged  intensity  are  developed.  Further,  in  Section  2.2,  the  theoretical 
advantage  of  controlling  intensity  is  investigated. 


2.1  Intensity  Expressions 

As  defined  in  Section  1.2.1,  the  total,  instantaneous  intensity  at  any  given 
point  is  the  sum  of  two  terms:  the  force  term  and  the  moment  term. 

The  force  tc>m  is  equal  to  Fv.  If  E  represents  the  Young’s  modulus,  and  I 
represents  the  moment  of  inertia  of  the  cross  section,  then  the  force  can  be  repre¬ 
sented,  in  terms  of  the  transverse  displacement  as 


F-£/. 


£i. 


The  velocity,  v,  can  be  expressed  as 


(2.1) 


(2.2) 


The  moment  term  is  equal  to  MQ.  The  moment  can  be  represented  by 
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dx^ 

The  angular  velocity,  Q,  can  be  expressed  as 


(2.3) 


Q=Ji. 

dtdx 


(2.4) 


The  total,  instantaneous  intensity,  therefore,  can  be  expressed  as 


dx^  dt  dx^  dtdx 


(2.5) 


It  follows  that  the  time-averaged  intensity  expression  given  in  Eq.  (1.2)  can  be 
represented  by 


<n>=i/?e 
‘  2 


\EI. 


d% 


dx' 


dt 


I-Irb 


Fl^ 

il' 

dx^ 

dtdx 

(2.6) 


2.2  Theoretical  Advantage  of  Controlling  Intensity 

The  following  developments  were  performed  by  Dr.  S.D.  Sommerfeldt  (1991) 
in  an  Acoustics  class  at  the  Pennsylvania  State  University. 

Consider  the  problem  of  controlling  the  vibration  in  an  infinite  beam,  as  shown 
in  Figure  2.1.  It  is  assumed  that  the  beam  is  excited  by  a  driving  force,  located 
at  x=0,  and  a  complex  control  force. 
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<2-7) 

located  at  x=L.  F„  is  taken  to  be  real  with  no  loss  in  generality.  The  region  of 
interest  is  taken  to  be  xiL.  The  two  solutions  considered  are  the  control  of  time- 
averaged  intensity  and  the  control  of  the  transverse  acceleration  at  a  point  x^. 


x-0  x-L 


Figure  2.1  Diagram  of  Infinite  Beam  with  Control. 


The  control  of  time-averaged  intensity  is  considered  first.  The  displacement 
of  the  beam  can  be  represented  by 


_ LJ-F  e  ->^*jF  e  ***-/  e  ->^-*KjF  e 

4EIk^'-  "  ■'  *  '  '  y 


i-^-F^e  -^*jFe  -^-Fe 


AEIk^ 


where  ;=v^,  and  k  is  the  flexural  wavenumber.  Internal  structural  losses  have  been 
ignored  in  this  formulation.  From  Eq.  (1.2),  the  time-averaged  intensity  at  any  point 
can  be  expressed  as 
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<n>=l/?e{Fv*}+l/?e{A/Q*}.  (2.9) 

'  2  2 

In  order  to  calculate  the  control  solution,  one  must  first  determine  the  intensity 
expression  for  Performing  the  necessary  partial  differentiation  results  in 

expressions  for  F,  v,  M,  and  Q: 


and 


(2.10) 


(2.11) 


(2.12) 


Substituting  these  expressions  into  Eq.  (2.9)  gives 


(2.13) 
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(2.14) 

In  order  to  determine  the  desired  control,  the  time-averaged  intensity  is  minimized 
with  respect  to  the  control  force.  Because  the  control  force  is  complex,  this  involves 
two  partial  differentiations: 


(2.15) 


(2.16) 


The  solution  of  the  two  equations  (2.15)  and  (2.16)  yields  the  requisite  control  force 
components: 


F^-F^osm 


(2.17) 


(2.18) 


Using  Euler’s  relation,  these  can  be  combined  to  yield 


(2.19) 


Substituting  this  into  Eq.  (2.14)  gives 
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{n>,=0;jciL 


(2.20) 


Further,  substituting  Eq.  (2.19)  into  Eq.  (2.8)  yields 


AEik^ 


(2.21) 


that  is,  decays  exponentially  to  0  for  x>L.  These  relations  are  true  regardless 
of  the  location  of  the  intensity  measurement  used  for  control. 

An  analogous  development  is  given  for  the  control  of  the  acceleration  at  a 
single  point.  Again  considering  the  infinite  beam  configuration  shown  in  Figure  2.1, 
the  acceleration  for  can  be  obtained  from  Eq.  (2.8)  as 


a{x,t)^Sl:^-Fe  -***jFe  -^-Fe  e 


This  acceleration  is  set  equal  to  0  at  some  location  Xo>L: 


Solving  for  F^  yields 


(2.22) 


(2.23) 


F»-_— if_F. 


(2.24) 


Substituting  this  into  the  expression  for  displacement  yields  an  expression  which,  in 
general,  is  not  equal  to  0  for  x^.  If,  however,  the  point,  x,,  is  far  removed  from  L, 
such  that  ibCo>>l  and  k(xQ-L)»l,  the  control  force  reduces  to 
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F  =-e  ->^F .  (2.25) 

C  M 

This  is  the  same  control  force  that  is  found  in  Eq.  (2.19)  for  the  case  of  controlling 
the  time-averaged  intensity  at  any  point  jn&L.  As  is  the  case  with  intensity  control, 
substituting  this  control  force  into  the  expression  for  displacement  at  results  in 
%{x,t)  decaying  exponentially  to  0  for  x>L. 

The  problem  with  controlling  acceleration  is  that  the  error  sensor  must  be 
located  far  (relative  to  a  wavelength)  from  the  source  or  any  boundaries  or 
discontinuities.  In  any  but  the  simplest  structures,  this  is  difficult  if  not  impossible. 

The  theoretical  advantage  of  controlling  intensity,  time-averaged  at  least,  is 
that  the  error  sensor  can  be  located  as  close  to  the  control  source  as  geometrical 
constraints  allow.  With  complicated  structures,  such  as  frames,  this  suggests  that  the 
vibration  could  be  confined  to  the  region  between  the  vibration  source  and  the  control 
source— entirely  within  the  first  spar. 


Chapter  3 

CALCULATION  OF  INTENSITY 

Many  different  ways  of  determining  structural  intensity  have  been  studied  in 
the  past  few  years.  Sensors  ranging  from  accelerometers,  strain  gauges,  and  force 
gauges  to  laser  velocimeters  and  vibrometers  have  been  employed  in  this  endeavor. 
In  this  study,  accelerometers  alone  are  considered  for  the  measurement  of  total, 
instantaneous  intensity.  This  does  not,  however,  determine  the  method  for  utilizing 
these  point  acceleration  measurements  to  calculate  the  structural  intensity. 

An  array  of  5  accelerometers,  positioned  on  the  top  of  the  beam  is  used.  The 
spacing  between  accelerometers  is  constant.  As  is  shown  in  Figure  3.1,  the 
accelerometers  are  arrayed  in  one  dimension,  and  the  accelerometer  positions  are 
numbered  sequentially.  The  acceleration  signal  at  position  i  at  time  t  is  referred  to 


BEAM 


Figure  3.1  Accelerometer  Array  Positioning. 
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In  Eq.  (2.5)  the  total,  instantaneous  intensity  is  expressed  as  a  function  of  the 
transverse  displacement.  If  a  time-harmonic  vibration,  with  frequency  (Oq,  is  assumed, 
then  the  displacement,  as  a  function  of  the  transverse  acceleration,  can  be  represented 
by 

(3.1) 

(Oo 

With  this  expression  for  the  total,  instantaneous  intensity  at  position  3  is  expressed 
as 


n3=£/ 


(  -1 


\/ 


<«>b 


2  dx^ 


-1  ^^3^ 


<4 


-El 


(  -1 


\/ 


tOo 


2  dx^ 


-1 


a>o 


(3.2) 


In  order  to  calculate  this  expression  for  structural  intensity,  the  various  time  and  space 
derivatives  or  the  acceleration  at  position  3  must  be  calculated. 

This  chapter  describes  the  various  calculation  tools  used  in  this  study  to 
determine  the  necessary  derivatives  and  the  resulting  structural  intensity.  The  chapter 
also  describes  some  additional  strategies  that  were  studied  and  discarded.  Section  3.1 
describes  the  technique  used  to  calculate  the  time  differentials  of  the  accelerations; 
Section  3.2  describes  the  technique  used  to  calculate  the  space  differentials  of  both 
the  accelerations  and  the  previously  calculated  time  differentials.  Section  3.2  also 
shows  the  conclusion  of  the  intensity  calculations.  Section  3.3  describes  the 
importance  and  the  method  of  calibrating  the  accelerometers  in  the  intensity  array. 
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Section  3.4  describes  several  other  techniques  studied  for  dealing  with  time 
differentials. 

3.1  Time  Differentials 

The  technique  used  to  calculate  the  time  differentials  of  the  accelerations 
(these  time  differentials  will  be  used  to  determine  the  velocity  and  angular  velocity 
terms)  is  a  backward  Hnite  difference  technique.  The  particular  finite  difference 
scheme  chosen  is  that  with  error  on  the  order  of  where  x  is  the  time  increment 
between  samples.  The  scheme  is  chosen  to  be  a  backward  finite  difference  scheme 
because  the  accelerations  are  known  only  up  to  and  including  the  time  at  which  the 
time  derivatives  are  desired.  The  scheme  is  chosen  to  be  that  with  an  error  on  the 
order  of  x^  because  the  error  at  this  step  must  be  minimized  to  keep  the  error  low  in 
the  ensuing  spatial  derivative  calculations. 

The  velocity  term,  v,  at  time  t  is  obtained  by  implementation  of  the  backward 
finite  difference  formula  on  acceleration  a,  at  time  t,  that  is  a^y. 


2x 


+0(x^. 


(3.3) 


In  order  to  determine  the  angular  acceleration  term,  Q,  at  time  r,  it  is  necessary 


to  calculate  the  time  differentials  at  positions  2  and  4 
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dt  2t 


'*’^2./ -2 


+0(t2) 


(3.4) 


and 


dt 


2t 


+0(t^). 


(3.5) 


3.2  Space  differentials 

Space  differentials  are  used  to  calculate  the  force,  moment,  and  angular 
velocity  terms.  The  technique  used  to  calculate  the  space  differentials  is  a  central 
finite  difference  scheme.  The  particular  scheme  has  error  on  the  order  of  A^,  where 
A  is  the  spacing  between  accelerometers.  This  scheme  was  chosen  because  Suen  [17] 
has  shown,  that  of  the  schemes  that  he  studied,  this  "second  order  scheme  ...  caused 
the  least  error  in  computing  the  structural  intensity." 

The  force  term,  F,  requires  the  calculation  of  the  third  spatial  derivative  of  the 
acceleration  at  point  3.  (Note:  The  lack  of  a  time  subscript  implies  the  term  is  at 
time  t.)  Using  a  central  difference  formula,  this  spatial  derivative  can  be  obtained  as 


d’fl-  -a,+2a,-2a.+a, 

_ 1=— 1 - 1 — L_l40(A^. 

dx^  2A’ 


(3.6) 


The  moment  term,  M,  requires  the  calculation  of  the  second  spatial  derivative 


of  the  acceleration  at  point  3,  which  can  be  obtained  as 
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The  angular  velocity  term,  Q,  requires  the  calculation  of  the  first  spatial 
derivative  of  the  time  differential  of  Oy  This  scheme  employs  the  time  derivatives 
of  flj.*  and  0^1  shown  previously  in  Eqs.  (3.4)  and  (3.5).  This  spatiai  derivative  is 
obtained  as 


-da^  da^ 


d  .  dt  dt 
drdx  “  Ia 


40  (A^. 


Assuming  the  above  partial  derivatives  have  been  calculated,  the  total, 
instantaneous  in.r‘'isuy  can  be  calculated.  Due  to  the  fact  that  the  operations  have 
y  performed  or  celerations,  and  not  displacements,  a  double  integration  must  be 
perfonped.  This  double  integration  is  achieved  by  multiplying  each  term  by  -l/(0o^ 
where  uig he  frequency  of  excitation.  (This  implies  that  this  method  of  calculating 
intensity  is  only  valid  for  single  frequencies,  with  the  frequency  of  excitation  known.) 
Thus  the  quantities  needed  to  determine  the  structural  intensity  can  be  summarized: 


El 


_d|3_  1  dflj 

a/  dt  ’ 


(3.10) 
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and 


A/3=-£/. 


f!i 

dx^ 


El 


Oil 


(3.11) 


Q  =^=-J_^. 

^  dtdx  ta! 

WQ 


(3.12) 


Because  F  and  v  are  multiplied  and  M  and  Q  are  multiplied,  the  resulting  factor  is 
1/oiQ*.  Realizing  that  both  F  and  M  have  factors  of  El,  the  resulting  expression  for 
total,  instantaneous  intensity  is  calculated  as 


^EI 

A 

«0 

dx^  dt  dx^  dtdx 

(3.13) 


3  J  Calibration  of  the  Accelerometers 

When  applying  spatial  Hnite  differences,  it  is  important  to  calibrate  the 
accelerometers  relative  to  a  reference  accelerometer.  The  accelerometers  do  not 
necessarily  have  to  be  calibrated  absolutely  because  the  adaptive  nature  of  the  control 
algorithm  will  correct  for  any  uniform  bias  across  the  array. 

Due  to  the  single  tequency  nature  of  the  excitation,  the  accelerometers  are 
calibrated  at  speciflc  frequencies  of  interest.  To  do  this,  two  accelerometers  (one 
reference  accelerometer  and  one  to  be  calibrated)  are  excited  with  an  impulse  and  the 
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response  of  each  recorded.  Then  a  ratio  of  the  magnitude  of  the  responses  is  taken 
at  various  frequencies.  For  example,  if  the  reference  accelerometer  is  referred  to  as 
R,  then  the  calibration  coefficient  of  accelerometer  A  at  300  Hz  is  equal  to  the 
magnitude  of  the  response  of  R  at  300  Hz  divided  by  the  magnitude  of  the  response 
of  A  at  300  Hz. 

3.4  Additional  Schemes  Studied 

This  study  investigated  other  strategies  for  dealing  with  the  time  dependence 
in  the  calculation  of  intensity.  These  dealt  with  the  problem  differently.  The  spatial 
finite  di^erence  scheme  was  used,  but  instead  of  taking  time  differentials  and  then 
double  integrating  with  it  was  attempted  to  perform  a  numerical  integration  of 
the  accelerations.  A  single  integration  of  Cj,  a^,  and  a^  would  be  necessary  to  get  the 
velocity  term  and  to  get  the  velocities  needed  to  calculate  the  angular  velocity  term. 
Double  integration  of  all  five  accelerations  would  be  needed  to  get  the  displacements 
needed  to  calculate  the  force  and  moment  terms.  Several  numerical  integration 
techniques  were  studied,  including  various  Adams-Moulton  schemes  and  schemes 
adapted  from  digital  filter  integrators  [20]. 

By  means  of  simulating  the  integration  techniques,  it  was  found  that  all  of 
these  additional  schemes  gave  intolerable  error.  The  most  common  type  of  error  was 
a  linearly  increasing  error.  This  error  quite  often  started  out  small  and  grew  slowly. 
However,  the  goal  is  to  run  the  control  indefinitely,  and  any  continuing  growth  in  the 
error  is  deemed  intolerable.  Another  consideration  with  these  numerical  integration 
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techniques  is  picking  proper  initial  conditions.  The  question  arises  as  to  what  a  valid 
assumption  for  the  initial  velocities  and  displacements  is.  This  question  was  not 
resolved  due  to  the  fact  that  the  techniques  were  unacceptable.  One  advantage  with 
differentiation  is  that  one  does  not  need  to  make  any  assumption  about  the  initial 
conditions,  although  differentiation  does  introduce  high-frequency  noise  into  the 


results. 


Chapter  4 

DESCRIPTION  OF  CONTROL  ALGORITHMS 


Having  decided  upon  a  means  of  calculating  the  total,  instantaneous  intensity, 
the  various  means  of  controlling  this  intensity  are  investigated.  The  algorithms  used 
for  the  control  system  are  the  projection  algorithm,  which  is  used  for  system 
identification,  and  the  filtered-x  least-mean-squares  algorithm,  which  is  used  to  update 
the  control  filters.  These  algorithms  are  described  in  this  chapter.  The  system 
identification  algorithm  is  described  first,  in  Section  4.1,  because  it  is  necessary  to 
have  these  filters  converge  to  their  stable  values  before  the  control  filters  can 
converge  to  their  stable  values.  Section  4.2  describes  the  filtered-x  least-mean-squares 
algorithm. 

4.1  System  Identification 

The  algorithm  used  for  system  identification  is  called  the  projection  algorithm 
or  the  normalized  least-mean-squares  (NLMS)  algorithm.  In  this  application  of  the 
algorithm,  there  is  assumed  to  be  one  error  sensor,  measuring  the  error  signal,  E(m), 
where  m  is  the  discrete  time  index.  The  system  identification  algorithm  attempts  to 
find  functional  transfer  functions  correlating  the  input  trace  and  the  control  signals  to 
this  error  signal.  The  error  signal  can  be  considered  as  consisting  of  two  parts:  the 
signal  which  is  present  in  the  absence  of  control  {le.  that  which  one  is  trying  to 
control)  and  that  which  arises  due  to  control  forces  driving  the  system.  In  the  early 
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literature  [11:  9-10]  the  element  of  the  error  signal  one  is  trying  to  control  is  referred 
to  as  the  "desired  signal,"  d(m);  apparently,  because  this  is  the  signal  one  "desires" 
to  match  with  the  adaptive  system.  If  one  assumes  that  N  control  actuators  are  used, 
the  element  of  the  error  signal  which  arises  due  to  control  forces  added  to  the  system 
can  be  approximated  by  passing  the  control  outputs,  y,(/n),  through  corresponding 
transfer  functions,  HJim).  The  transfer  function  from  the  /i*  control  actuator  to  the 
error  signal  is  represented  by  a  linear  Finite  Impulse  Response  (FIR)  frlter  vector  with 
K  coeffrcients,  or  taps, 

(4-1) 

The  time  train  of  y,  necessary  to  implement  this  FIR  filter  also  has  K  terms  and  is 
represented  by 

[y^(m)y^(/n-l)-y,(m-#:+l)].  ('♦•2) 

The  element  of  the  error  signal  that  arises  due  to  the  n*  control  actuator  can  be 
represented  by  a  convolution  sum 

In  vector  notation,  this  is  expressed  as 

For  notational  convenience,  a  composite  control  path  transfer  function  vector  can  be 
represented  as  a  concatenation  of  the  N  transfer  function  vectors. 
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Likewise,  a  composite  control  output  vector  can  be  formed  as  a  concatenation  of  the 

N  time  trains, 

Y'(m)m  [K,''(m)|r/(m)|-|y/(m)).  (‘‘•6) 

Using  this  notation,  the  total  error  signal  can  be  expressed  as 

e(m)^d(m)*Y  \m)H(m).  (4- ”7) 

However,  in  implementing  the  projection  algorithm,  a  necessary  condition  to  achieve 
convergence  is  that  the  entire  error  signal  be  approximated.  A  necessary  condition 
for  the  control  algorithm  to  provide  attenuation  is  that  the  desired  signal,  d(m),  be 
correlated  with  the  input  trace.  This  relationship  can  be  modelled  through  the 
implementation  of  an  additional  transfer  function  Alter 

If  the  x(m-k)  vector  is  concatenated  with  l^(m),  one  obtains  the  data  vector. 

Likewise,  if  the  vector  is  concatenated  with  tf{m\  one  obtains  the 

coefficient  vector. 
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e'W-  Vw-,,]- 

Now  the  total  error  signal  used  by  the  system  identiflcation  algorithm  can  be 
expressed  as 

eim)‘ Q^imyp(m)  (^-H) 

[21:  70-75].  One  can  calculate  y,(m),  and  one  can  measure  E(m)  and  x(m).  With 
this  knowledge,  the  Projection  (or  NLMS)  Algorithm  can  be  employed  to  converge 
upon  values  for  6(/r).  If  the  estimated  value  of  6(m)  is  represented  by  ^m),  then 
this  estimate  is  updated  according  to  the  Projection  Algorithm. 

^in+l)=  ^/n)+ - - [e(m)- 6^(m)^m)l  (^>12) 

b^(myp(my  ^ 

where  is  the  updated  transfer  function  estimate,  a  is  a  convergence 

parameter-chosen  to  be  0<a<2  for  stability-and  b  is  a  small,  positive  constant  used 
to  ensure  that  there  is  no  division  by  zero  [22:  50-54]. 


4.2  Control 

The  algorithm  used  for  the  control  Hlter  update  is  called  the  filtered-x  least- 
mean-squares  (filtered-x)  algorithm;  it  is  similar  to  one  developed  by  Widrow  and 
Steams  [11:  288-294].  This  algorithm  is  used  to  determine  the  N  control  signals, 
yj(m),  which  minimize  a  mean-squared-error  signal.  The  control  signals  are  obtained 
by  passing  the  system  input,  x(m),  through  the  control  filters,  W,(m).  These  control 
filters  are  each  represented  by  a  FIR  filter  vector  with  I  taps 
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The  time  train  of  the  input  signal,  jc(m),  necessary  to  implement  this  FIR  filter  also 
has  /  terms  and  is  represented  by 


X\m)m  [jc(m)jr(m-l)-x(w-/-^l)]. 

A  convolution  sum  of  these  two  signals  gives  the  n*  control  output 


(4.14) 


Ml 

where  wj^m)  are  the  n*  control  filter  coefficients  at  time  m,  y,(m)  is  the  n*  control 
filter  output  at  time  m,  and  x(m-i)  is  tl«  input  to  the  control  filters  at  time  m-i.  A 
composite  control  filter  vector  can  be  formed  as  the  concatenation  of  the  N  control 
filter  vectors: 


The  optimal  solution  for  the  control  filter  vector  is  taken  to  be  the  solution 
which  minimizes  the  mean-squared-error  signal.  The  filtered-x  algorithm  is  a  steepest 
descent  algorithm  which  is  based  on  the  concept  of  updating  the  control  filter 
coeffiewnts  according  to  the  negative  of  tlw  gradient  with  respect  to  the  filter  coeffi¬ 
cients  of  the  mean-squared  error,  so  that  the  control  filters  converge  to  the  optimal 
solution  in  an  iterative  manner.  The  mean-squared  error  is  approximated  by  the 
instantaneous  squared  error,  e^(m).  This  error  signal,  e(m),  is  not  necessarily  the  same 
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as  the  error  signal  used  in  the  system  identification  algorithm,  t(m).  The  update 
equation  is  expressed  as 

(4.17) 

The  so-called  convergence  parameter,  Po*  is  a  parameter  chosen  to  ensure  convergence 
and  stability.  The  squared  error,  e^(m),  in  least-mean-squares-based  algorithms,  is 
desired  to  be  quadratic  in  W  [11:  19-21].  In  the  acceleration  control  application, 
e{m)  is  equal,  therefore,  to  the  acceleration.  However,  in  the  intensity  control  case, 
the  intensity  is  quadratic  in  terms  of  W.  This  is  because  F,  v,  M,  and  Q  are  each 
linear  in  W.  Therefore,  the  force  term,  Fv,  and  the  moment  term,  AfQ,  and, 
subsequently,  the  intensity  are  each  quadratic  in  W.  In  order  to  have  e^(m)  quadratic 
in  W,  e(m)  is  set  equal  to  the  square  root  of  intensity.  If  11  represents  the  instanta¬ 
neous  intensity,  then  the  update  equation  can  be  expressed  as 

W(m+l)-W(m)-|i„V^.  (4.18) 

The  explicit  forms  of  the  gradients  are  different  for  each  of  the  various  control 
approaches  used  in  this  thesis.  Therefore,  the  developments  of  the  gradients  follow 
in  the  appropriate  sections. 

In  all  control  approaches  used,  it  is  necessary  to  determine  an  appropriate 
value  for  the  convergence  parameter,  pg.  Sommerfeldt  [21:  84-89]  has  shown  that, 
for  the  filtered-x  algorithm,  a  bound  on  pg-necessaiy  for  convergence-is  that 
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0<2vi,<^,  (4.19) 

where  is  the  largest  eigenvalue  of  the  auto-correlation  matrix, 
where  R(m)  is  the  input  signal,  X(m),  filtered  by  the  composite  transfer  function  filter, 
H{m).  If  the  system  has  1  error  sensor  and  N  control  actuators,  the  control  niters 
have  /  taps,  and  the  transfer  function  filters  have  K  taps,  then  an  approximation  for 
this  expression  can  be  obtained  using  the  result 

where  RJO)  is  the  average  power  of  the  input  signal,  and  is  the  maximum  value 
of  the  control  path  transfer  functions*  coefficients.  Thus  a  conservative  range  for  Po 
can  be  given  by 


0<2po<-_ - 1 - 

hLUf*I*K^]RJO) 


(4.21) 


Adaptive  implementation  of  this  can  be  facilitated  by  calculating  a  value  of  the 
average  input  power,  by  monitoring  the  transfer  function  coefficients,  and  by, 
subsequently,  updating  the  value  of  Po  at  each  time  step. 


4.2.1  Control  of  Acceleration 

When  controlling  acceleration,  the  filtered-x  algorithm  can  be  applied  in  a 
straightforward  manner.  The  error  signal  is  the  acceleration  one  seeks  to  control. 
This  acceleration  is  chosen  to  be  that  at  position  3  in  Figure  3.1,  Oy  In  this  thesis. 
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one  control  actuator  is  used  to  control  the  acceleration.  There  is,  therefore,  a  single 
control  filter  and  a  single  control  path  transfer  function  filter.  Due  to  the  single 
actuator,  the  subscript  n  is  dropped.  The  error  signal,  e(m),  can  be  represented  by 

e(m)=d(m)*H\m)Y(m).  (^-22) 

In  summation  notation  this  is  expressed  as 

K-l 

e(m)=d(m)+^  hjjn)y{m-k).  (^-23) 

Using  the  FIR  control  filter  definition  this  becomes 

kMi  iaO 

The  standard  development  utilizes  the  assumption  that  w^m)  is  time  invariant,  which 
allows  the  order  of  the  summations  to  be  interchanged.  If  this  interchange  is  made, 
the  error  signal  can  be  expressed  as 

e{m)^m)*Y,  h^{m)x{m-i-k).  (^-25) 

If  r(m)  is  defined  to  be 

(^-26) 

then  the  error  signal  can  be  expressed  as 
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i-i 

g(m)=fi(m)-»y^  w.  r(m-i).  ('*•27) 

i«0 

It  can  be  seen  that  the  signal  r(m)  corresponds  to  the  input  signal,  x(m),  being  Altered 
by  the  transfer  function  H,  which  leads  to  the  name  filtered-x.  Furthermore,  if  the 
vector  R(m)  is  defined  as 

R  \m)m  [r(m)r(m-l>-r(m-/-^l)], 
then  the  error  signal  can  be  rewritten  in  vector  notation  as 

e(m)=d(m)-^W^R. 

From  this  expression  the  square  of  the  error  is  found  to  be 


(4.28) 


(4.29) 


e\m)Md\m)*2d(m)R  ^(m)W(m)*W  ^(m)R(m)R  \m)W{m). 
The  gradient  of  this  result  can  be  expressed  as 


(4.30) 


V^*-2R(m)d(m).^2R(in)«'’(m)W(m)  (4-31) 

or 


V^*«2R(m)[d(in)>»«  '■(in)W(m)]«2R(m)e(w).  (4*32) 

From  this  result,  the  update  equation  for  the  control  filter  coefficients  can  be 
expressed  as 

W(m*l)-HXifi)-pe(m)J?(m),  (^-33) 


where 
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fi  =  2|io 

is  a  parameter  chosen  to  ensure  convergence  and  stability,  and  e{m)  is  the  measured 
acceleration  signal.  One  should  note  that  the  error  signal  used  for  system  identifica¬ 
tion,  E(m),  is  also  this  acceleration;  that  is, 

e(m)  =«(/«)  (4-34) 

This  is  because  the  filtered  quantity  that  the  system  identification  algorithm  is  trying 
to  model,  R,  corresponds  to  this  acceleration. 

4.2.2  Control  of  Intensity 

The  intensity  control  application  of  the  filtered-x  algorithm  is  not  as 
straightforward  as  the  acceleration  control.  This  is  because  the  total,  instantaneous 
intensity  is  a  quadratic  quantity  in  W.  Thus,  to  take  the  gradient  of  the  squared 
intensity  would  not  be  a  valid  application  of  the  theory.  Instead,  the  gradient  with 
respect  to  the  control  filter  coefficients  is  taken  of  the  intensity: 

W(m+l)=W(w)-p^V,^n  (4.35) 

or 

(4.36) 

Off 

To  obtain  an  expression  for  the  gradient,  the  intensity  can  be  represented  by 
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n=Fv+A/Q. 

One  should  note  that 


and  M^-Ei£!L,  (4.38) 

dx^  dx^ 

where  |  is  the  transverse  displacement.  This  allows  the  intensity  to  be  expressed  as 


n=£/ 

_L_  (vdt*v 

-El 

dxO 

dxJ 

Use  has  also  been  made  of  the  fact  that 


(4.39) 


(4.40) 

dx 

Further,  if  v  and  Q  are  separated  into  a  desired  component  (the  response  with  no 
control)  and  a  component  due  to  the  control,  they  can  be  represented  by 


(4.41) 


By  interchanging  the  order  of  summation,  as  described  in  Section  4.2.1,  v  and  Q  can 
also  be  expressed  as 


v«Vg+W 


(4.42) 


Substituting  these  expressions  into  Eq.  (4.39),  the  intensity  can  be  expressed  as 


V^=£/JL  [Rdt*v*F*R-ElA.  [RJt*Q^M*R^.  (4.45) 

dx^J  '  dxJ  °  “ 

This  expression  is  exact,  but  it  is  impractical  for  implementation  on  a  real-time  DSP 
board.  To  achieve  a  practical  control  scheme  for  implementation,  two  sets  of 
heuristic  assumptions  are  used  to  obtain  approximations  of  this  gradient.  For  the  first 
approximation,  the  assumptions  are  made  that  the  third  partial  derivative  with  respect 
to  X  yields  a  factor  of  that  the  first  partial  derivative  with  respect  to  x  yields  a 
foctor  of  'jkj  and  that  the  integral  with  respect  to  time  yields  a  factor  of  (j  is 
again  k  is  again  the  flexural  wavenumber,  and  (Uq  is  the  driving  frequency.)  It 
should  be  noted  that  using  these  assumptions  corresponds  to  assuming  time-harmonic 
excitation  and  that  implementing  the  resulting  algorithm  requires  a  knowledge  of  the 
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driving  frequency.  Using  these  assumptions  in  equation  (4.45)  allows  the  gradient  to 
be  approximated  by 

(4.46) 

“o  "o 

One  of  the  possible  limitations  is  that  the  signs  of  the  various  terms  depend  on  the 
type  of  wave  which  is  considered.  The  signs  chosen  here  correspond  to  a  propagating 
wave  in  the  -i-x  direction.  Obviously,  if  another  type  of  wave  dominates  the  vibration 
field,  some  of  the  phases  in  the  gradient  will  be  wrong,  and  that  may  affect  the 
performance  of  the  control  system.  Because  Hq  and  J?„  correspond  to  an  angular 
velocity  and  a  transverse  velocity,  respectively,  they  are  related  according  to: 

(4.47) 
dJC 

If  the  derivative  is  again  replaced  by  -jk,  R^  can  be  approximated  as 

Ra—JkR, 

These  results  allow  the  gradient  to  be  approximated  as 


EIk\  _  .Elk^  ^ 
*v*F-j - *Q  ~jkM 


a)„ 


(4.49) 


This  expression  cannot  be  implemented  using  the  real  values  obtained  from 
accelerometers.  If  the  various  signals  are  assumed  to  be  in  phase~an  arguable 
assumption  at  best,  and  if  the  -jk  is  replaced  by  k,  then  the  gradient  can  be 
approximated  by 
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.  (^-50) 

[%  Wo 

If  this  expression  is  used  in  the  filtered-x  algorithm,  the  update  equation  can  be 
expressed  as 

W(m+l)=W(m)-pe(m)R,,  (4.51) 

where 

e(m)=fiL*v+F+£lfL*Q4*W  (4.52) 

"o  Wo 

is  the  effective  error  signal  for  the  control  algorithm.  It  should  be  noted  that  this  is 
not  the  effective  error  signal  for  the  system  identification  algorithm.  In  this  case  the 
error  signal  for  the  system  identification  algorithm  is  the  velocity  term,  v.  This  is 
because  the  filtered  quantity  that  this  algorithm  is  trying  to  model  is  the  filtered 
quantity  correlating  the  control  output  to  the  velocity.  So  for  this  representation  of 
the  gradient  the  error  signal  used  in  the  system  identification  is  expressed  as 

e(m)»v.  (4*53) 

A  simplification  of  Eq.  (4.49)  can  be  obtained  by  again  assuming  a  time- 
harmonic,  propagating  wave  in  the  -t-x  direction.  In  this  case, 

P  Elk^  .  w  Elk^ 

F*Et — 4* - V  and  M^-EI — 2.* - Q. 

dx^  <Op  dx^  «o 

(4.54) 

Using  these  results  in  Eq.  (4.49)  allows  the  gradient  to  be  approximated  by 
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(4.55) 

"o 

The  update  expression  can,  therefore,  be  expressed  as 

W(/n+l)=W(m)-pe(m)R,,  (4.56) 

where 

e(«)=ifcv+fl  (4-57) 

is  the  effective  error  signal  for  the  control  algorithm.  Again,  it  should  be  noted  that 
this  is  not  the  effective  error  signal  for  the  system  identiHcation  algorithm.  The 
system  identification  error  signal  is  the  velocity;  that  is, 

«(m)=v.  (4.58) 

This  is  because  the  algorithm  is  trying  to  model  the  filtered  quantity, 

The  terms  used  in  the  calculation  of  the  two  intensity  error  fimctions,  Eqs. 
(4.52)  and  (4.57),  are  the  same  terms  used  in  the  calculation  of  the  intensity.  For 
example,  the  velocity  term,  v,  used  in  determining  both  error  function  signals  is 
calculated  as  in  Eq.  (3.10),  described  in  Sections  3.1  and  3.2.  It  should  be  noted  that 
the  intensity  signal  is  not  directly  used  in  the  implementation  of  these  control 
algorithms,  but  this  signal  is  calculated  as  a  basis  of  evaluating  the  effectiveness  of 
the  control. 

In  addition  to  several  possible  error  signals,  three  different  actuator 
configurations  are  considered.  These  consist  of  controlling  the  beam  with  a  control 


37 


force,  a  control  moment,  or  both  a  force  and  a  moment.  Because  there  are  two 
components  to  intensity— force  and  moment— the  force-moment  configuration  for 
control  is  believed  to  be  the  most  promising.  The  implementation  of  each  of  these 
three  strategies  is  discussed  briefly. 


i 


4.2.2. 1  Control  of  Intensity  with  a  Force  Output 

The  force  output  is  calculated  with  one  control  filter.  There  is,  therefore,  a 
single  control  filter  and  a  single  control  path  transfer  function  filter.  The  force  output 
is  sent  to  two  control  shakers.  This  is  done  in  such  a  way  that,  if  the  output  to 
control  shaker  1  is  represented  by  the  output  to  control  shaker  2  is  represented  by 
and  the  force  output  is  represented  by  yi,  then  y^  and  y^  are  expressed  as 


It  is  assumed  that  the  control  shakers  are  positioned  close  to  each  other-relative  to 
a  wavelength— so  that  the  effect  on  the  beam  approximates  a  force  applied  at  a  point 
halfway  between  the  shakers.  Although  there  are  two  shakers,  there  is  considered  to 
be  only  one  control  actuator;  this  is  done  to  make  the  development  analogous  to  that 
of  the  moment  control. 


4.2.2.2  Control  of  Intensity  with  a  Moment  Output 

The  moment  output  is  calculated  with  one  control  Alter.  There  is,  therefore, 
a  single  control  filter  and  a  single  control  path  transfer  function  Alter.  The  moment 
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output  is  sent  to  one  control  shaker,  and  the  negative  of  the  moment  output  is  sent  to 
the  other  control  shaker.  That  is,  if  the  output  to  control  shaker  1  is  represented  by 
y^i,  the  output  to  control  shaker  2  is  represented  by  y^,  and  the  moment  output  is 
represented  by  yj,  then  and  y^  are  expressed  as 


(4.60) 


It  is  assumed  that  the  control  shakers  are  positioned  close  to  each  other— relative  to 
a  wavelength-so  that  the  effect  on  the  beam  approximates  a  moment,  or  couple, 
applied  at  a  point  halfway  between  the  shakers.  Again,  although  there  are  two 
shakers,  there  is  only  one  control  actuator. 


4.2.2.3  Control  of  Intensity  with  Force  and  Moment  Output 

The  force-moment  actuator  configuration  for  control  is  a  superposition  of  the 
force  configuration  and  the  moment  configuration.  The  force  output  is  calculated  with 
one  control  filter,  and  the  moment  output  is  calculated  with  another  control  filter. 
There  are,  therefore,  two  control  filters  and  two  control  path  transfer  function  filters. 
The  force  output  is  sent  to  both  shakers.  The  moment  output  is  sent  to  oim  shaker 
and  the  negative  of  the  moment  ouQ)ut  is  sent  to  the  other  shaker.  That  is,  if  the 
output  to  control  shaker  1  is  represented  by  y^,  the  output  to  control  shaker  2  is 
represented  by  y^j,  the  force  output  is  re{»esented  by  y„  and  the  moment  output  is 
represented  by  yj,  then  y^  and  y^  ate  expressed  as 


“>'»("») 
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(4.61) 


4.23  Summary  of  Control  Update  Equations 

The  update  equations  to  be  used  in  this  thesis  are  common  to  all  strategies  for 
control.  The  system  identiflcation  filter  update  equation  is  expressed  as 

^i»i+l)«  6(m)+ - ^5*^^ - [  e(m)-  ©^(m)<I>(m)l  (4.62) 

6'*^(m)^m) 

and  the  control  filter  update  equation  is  expressed  as 

W(m+l)*lV(m)-pe(m)Jl(m).  (4.63) 

The  variations  in  control  strategies  result  fiom  the  definition  of  the  error  signals  and 
from  the  number  of  control  filters  employed  in  the  strategy. 

In  the  acceleration  control  strategy,  the  error  signals  are  defined  as 


e(m)^(/n)-a,(m),  (4.64) 

and  the  vectors  in  the  update  equations  ate  defined  as 

[y(m)-y(m-/ir+l)x(m)-x(m-/C+l)],  (4.66) 

W^(m)-[w,(m>..>v^,.,,(m)],  (4.67) 


and 
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JI  *■(«)=  [r(m)-r(m-/+l)].  (4.68) 

In  the  various  intensity  strategies,  the  control  update  error  signals  are 
expressed  as 

e(m)«f:l!L*v+F*Zl!L*Q*kM  (4.69) 

"o  "o 

or 

c(m)=ifev+Q,  (^-70) 

and  the  system  identiflcation  update  error  signal  for  both  instances  is  expressed  as 

e(m)*v.  (4*71) 

For  the  force  actuator  configuration  and  for  the  moment  actuator  configuration,  the 
update  equation  vectors  have  the  same  form,  and  they  are  expressed  as 


6^(m)«  [ *„(m)~.  A^^.^(m)  hj,my 

(4.72) 

^(^)*  Jc(m)-x(m-#r+l)], 

(4.73) 

W  '■(m)»  [  »Vj(m>-  Wp., )(«)]. 

(4.74) 

and 

R  ^(m)»  [r(iii)~’r(m -/+!)]. 

(4.75) 

For  the  force-moment  actuator  configuration,  the  update  equation 

vectors  are 

expressed  as 
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(4.76) 

<I>^(/n)=  y^my-y^m  -^T+l)  x(my-x(m  -^+1)] , 

(4.77) 

W  ^(m)=  [  wjmy-  w,(,.,/m)  wjmy-  Wj^^.i/w)],  (4.78) 

and 

X  ^(m)=  [rj(m>-r,(m-/+l)r2(m>.  r2(m-/+l)].  (4-79) 


Chapter  5 

EXPERIMENTAL  METHOD 


The  experimental  method  is  described  in  this  chapter.  The  chapter  is  divided 
into  two  parts:  Section  5.1  discusses  the  experimental  setup  and  Section  5.2  discusses 
the  experimental  procedure. 

5.1  The  System 

The  system  is  divided  into  two  parts:  the  system  being  controlled  and  the 
controlling  system.  Section  5.1.1  describes  the  structural  system  being  controlled, 
including  the  beam,  the  terminations,  and  the  support  system.  Section  5.1.2  describes 
the  control  system  hardware. 

5.1.1  Beam  Description 

The  structure  to  be  controlled  in  this  study  consists  of  a  thin-walled,  square, 
hollow  steel  beam;  one  end  of  which  is  embedded  in  a  box  of  sand.  The  other  end 
is  supported  with  a  compliant  rubber  band,  which  is  festened  to  an  overtiead  support 
structure.  This  rubber  band  is  used  to  support  the  static  load,  and  it  is  assumed  that 
the  rubber  band  does  not  significantly  alter  the  dynamic  loading  on  the  beam. 

The  steel  beam  has  been  found  to  have  the  following  ^ometrical  and 
mechanical  properties:  density  (p)  =  7.70x10’  kg/m’;  cross  sectional  area  (A)  = 
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74.0x10'®  m*;  mcxlulus  of  elasticity  (£)  =  187x10’  Pa;  and  moment  of  inertia  of  the 
cross  section  (7)  =  2.71x10’ m^. 

One  end  of  the  beam  is  embedded  in  a  box  of  sand  in  order  to  introduce 
absorption  at  the  boundary  and,  hence,  reduce  the  magnitude  of  the  reflection  coeffi¬ 
cient  from  the  value  of  1.0  associated  with  a  purely  free  termination.  A  lossless  beam 
and  a  reflection  coefficient  of  1.0  at  the  boundary  results  in  a  standing  wave  field 
only  and,  therefore,  only  reactive  (or  imaginary)  intensity  is  generated.  Reducing  the 
reflection  coefficient  at  the  boundary  reduces  the  reflected  wave  and,  therefore,  allows 
the  propagating  wave  component  to  dominate  the  wave  field.  This  results  in  an  active 
(or  real)  intensity  component.  If  the  termination  is  anechoic,  corresponding  to  the 
reflection  coefficient  having  a  value  of  0.0,  then  the  intensity  is  entirely  active. 
Previous  studies  on  beams  with  the  same  properties  and  similar  terminations  have 
shown  that  the  reflection  coefficient  is  approximately  0.254+j0.044  at  5000  Hz  and 
0.536+j0.058  at  500  Hz  [17,23].  Therefore,  for  the  frequency  range  used  for  this 
study  (140  Hz  to  520  Hz),  a  combination  of  active  and  reactive  intensity  is  expected. 

The  beam  is  excited  at  its  free  end  with  a  Ling  Dynamic  Systems  (LDS) 
model  V203  magnetic  driven  vibrator.  A  sinusoidal  signal  is  generated  by  a 
Wavetech  model  75  arbitrary  waveform  generator.  This  signal  is  passed  through  a 
BGW  Systems  model  200  professional  power  amplifier.  The  signal  is  then  sent  to  the 
LDS  V203  shaker,  which  excites  the  beam  with  a  sinusoidal  vibration. 

The  LDS  V203  shaker  is  connected  to  the  beam  by  a  stinger  to  insure  the 
delivery  of  a  moment-free  transverse  force  to  the  beam.  The  stinger  is  relatively 
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flexible  in  its  transverse  direction,  in  order  to  minimize  the  transfer  of  any  moment 
to  the  beam  due  to  misalignment.  This  stinger  is  constructed  with  a  short  piece  of 
20-gage  steel  wire  soldered  into  holes  drilled  in  two  bolts.  The  bolts  are  the  same 
size  as  the  tapped  mounting  hole  in  the  shaker.  One  bolt  is  threaded  into  the 
mounting  hole  in  the  shaker,  and  the  other  bolt  is  glued  with  epoxy  to  the  beam. 

5.1.2  Control  System 

The  control  system  consists  of  three  parts:  the  input  sensors,  the  processor, 
and  the  output  actuators.  A  schematic  diagram  of  the  control  system  is  foimd  in 
Figure  5.1.  This  figure  shows  the  various  parts  of  the  control  system  and  the 
connection  paths.  Section  5.1.2.1  discusses  the  input  sensor  hardware;  Section  S.1.2.2 
discusses  the  processor  hardware;  and  Section  S.1.2.3  discusses  the  output  actuator 
hardware. 

5.1.2.1  The  Inputs 

The  six  input  signals  to  the  control  system  originate  from  5  accelerometers  and 
the  signal  from  the  signal  generator  which  drives  the  excitation  shaker.  The  signals 
from  the  accelerometers  are  used  to  calculate  the  error  signals  (the  error  signal  in 
acceleration  control  is  the  third  acceleration  signal).  The  signal  from  the  signal 
generator  is  used  as  the  control  filter  input,  x(m). 

The  five  accelerometers  used  to  calculate  the  error  signals  are  PCB  Piezo- 
tronics,  Inc.  (PCB)  model  336A  accelerometers.  The  serial  numbers  of  the 
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Figure  5.1  Schematic  Diagram  of  Control  System, 
accelerometers,  in  order  of  position  from  the  free  end,  are  #1136,  #1156,  #1140, 
#1141,  and  #1142.  The  response  calibration  of  each  accelerometer,  as  described  in 
section  3.3,  is  performed  relative  to  accelerometer  #1136.  The  complete  table  of 
calibration  constants  calculated  is  given  in  the  Appendix.  The  accelerometer  signals 
pass  through  a  conditioning  amplifier,  which  is  a  PCB  483B07  12-channel  amplifier. 
This  amplifier  is  chosen  for  compatibility  with  the  PCB  accelerometers.  The  control 
filter  input  signal  is  also  passed  through  the  conditioning  amplifier. 

After  being  amplified,  the  six  input  signals  are  sent  to  a  Spectrum  Signal 
Processing  Inc.  32  channel  analog  infnit  board.  This  board  consists  of  a  12  bit  A/D 
converter,  and  it  allows  a  4  channel  sample  and  hold  operation  [24].  The  input 
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board  has  a  dynamic  range  of  ±  2.5  V.  It  is  designed  to  interface  with  Spectrum’s 
line  of  Digital  Signal  Processing  (DSP)  boards. 

5.1. 2.2  The  Processor 

The  processing  loop  includes  the  input  board,  the  DSP  board,  the  host 
computer,  and  the  output  board.  The  input  board  is  the  Spectrum  32  channel  analog 
input  board  described  previously.  The  DSP  board  is  a  Spectrum  Signal  Processing 
Inc.  DSP  96002  board,  which  is  built  around  a  Motorola  XSP96002RC33  chip.  The 
processor  is  a  96-bit  floating  point  processor.  The  DSP  96002  communicates  with  the 
input  and  output  boards  over  a  16  bit  data  bus.  The  output  board  is  a  Spectrum 
Signal  Processing  Inc.  16  Channel  analog  output  card. 

The  Spectrum  boards-input,  output,  and  DSP-all  requite  a  host  PC.  The  host 
PC  used  is  a  Hewlett  Packard  Vectra  QS/20,  386-ba5ed  machine.  The  host  PC 
provides  power  for  the  boards,  as  well  as  a  platform  for  generating  the  code  to  be 
downloaded  to  the  DSP  board.  The  DSP  board  data  and  memory  registers  are  also 
monitored  via  the  host  PC. 

5.1.2.3  The  Outputs 

The  outputs  consist  of  the  necessary  control  output  signals.  For  the  case  of 
acceleration  control  there  is  only  one  such  output  signal.  For  the  case  of  intensity 
control  (force  only,  moment  only,  or  force  and  moment)  there  are  two  such  output 
signals. 


47 


These  output  signal(s)  are  sent  through  a  Spectrum  Signal  Processing  Inc.  16 
channel  analog  output  card.  The  output  board  is  a  12-bit  D/A  converter  and  allows 
16  channels  of  simultaneous  output  through  the  use  of  a  double  buffer  output  system. 
If  the  double  buffer  is  used,  the  signals  are  sent  to  a  buffer  and  at  the  proper  time  all 
of  the  signals  are  released  simultaneously.  The  output  board  has  a  dynamic  range  of 
approximately  ±  8.19  V  [25]. 

The  output  signal(s)  are  ampliHed  by  means  of  a  BGW  Systems  model  200 
professional  power  amplifier. 

After  amplification,  the  signals)  are  used  to  drive  the  control  shaker(s).  The 
control  shakers  are  LDS  model  V102  magnetic  driven  vibrators.  The  V102  shakers 
employ  stingers,  whose  construction  is  similar  to  that  of  the  V203’s  stinger.  The  only 
difference  is  that  the  bolt  size  is  different,  due  to  the  fact  that  the  tapped  hole  in  the 
V102  shaker  has  a  different  size.  Again,  the  bolt  on  one  end  of  the  stinger  is 
mounted  to  the  shaker,  and  the  bolt  on  the  other  end  of  the  stinger  is  glued  with 
epoxy  to  the  beam. 

5.2  The  Experimental  Procedure 

This  section  describes  the  experimental  procedure.  The  measurement  of  the 
structural  intensity  is  conducted  at  one  frequency,  300  Hz.  This  frequency  is  chosen 
because  it  is  thought  that  the  spacing  and  timing  regarding  the  finite  difference 
techniques  are  optimal  at  300  Hz.  The  adaptive  vibration  control  investigation  uses 
frequencies  grouped  around  300  Hz.  The  adaptive  vibration  control  study  is 
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conducted  at  six  frequencies,  three  of  which  are  found  to  be  resonances  of  the  beam 
system:  140  Hz,  296  Hz,  and  488  Hz.  Three  additional,  off-resonance  frequencies 
are  also  studied:  200  Hz,  340  Hz,  and  520  Hz.  Section  5.2.1  describes  the  structural 
intensity  measurements,  and  Section  5.2.2  describes  the  adaptive  vibration  control 
implementation. 

All  signal  analysis  is  done  using  a  HP  3566A  signal  analyzer,  which  is 
referred  to  as  the  analyzer.  The  HP  3566A  is  a  frequency  analyzer  that  needs  a  host 
pc.  The  host  pc  is  the  HP  Vectia  QS/20  described  previously.  This  analyzer  rises  the 
MS-DOS©  operating  system  and  a  Microsoft©  Windows-based  user  interface. 

5.2.1  Intensity  measurement 

This  section  describes  the  study  of  the  structural  intensity  measurement.  An 
attempt  is  made  to  coirelate  frequency  spectra  of  the  calculated  output  and  input 
structural  intensities,  in  order  to  verify  the  algorithm  developed  in  Chapter  3.  The 
output  intensity  refers  to  the  total,  instantaneous  intensity  measured  at  the  accelerome¬ 
ter  array  and  calculated  using  this  algorithm.  The  input  intensity  is  a  measure  of  the 
intensity  input  to  the  beam.  It  is  measured  at  the  excitation  shaker.  The  assumption 
is  made  that  at  the  excitation  shaker  the  intensity  is  dominated  by  the  force  term.  For 
this  reason,  the  force  term,  Fv,  of  the  input  intensity  is  measured,  using  a  force  gauge 
and  an  accelerometer.  The  experimental  configuration  is  as  shown  in  Figure  5.2. 

The  input  intensity  is  calculated  by  multiplying  the  force  at  the  input  by  the 
velocity  at  the  input.  The  force  at  the  input  is  found  using  a  PCB  208B  force  trans- 
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Figure  5.2  Experimental  Setup  for  Intensity  Measurement. 

ducer,  positioned  between  the  shaker  and  its  stinger.  This  signal  is  conditioned  using 
the  PCB  483B07  amplifier  described  previously.  The  signal  is  then  sent  to  the 
analyzer.  The  velocity  at  the  input  is  found  by  integrating  an  acceleration  signal  at 
the  input.  This  integration  is  performed  in  the  firequency  domain,  through  multiplying 
by  I/70).  The  acceleration  signal  is  found  using  a  PCB  303A13  accelerometer.  The 
acceleration  signal  is  conditioned  through  the  PCB  483B07  amplifier  described 
previously.  This  acceleration  signal  is  then  sent  to  the  analyzer. 

The  input  intensity  ideally  would  be  calculated  by  first  transforming  the 
acceleration  into  the  firequency  domain.  Integration  of  this  acceleration  signal  would 
then  be  performed  in  the  frequency  domain,  by  means  of  multiplying  by  l/Jfo.  Then, 
this  velocity  signal  would  be  transformed  to  the  time  domain.  This  time  sequence  of 
the  velocity  signal  would  be  multiplied  by  the  time  sequence  of  the  force  in  order  to 
obtain  a  time  sequence  of  the  input  intensity. 
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However,  the  vibration  at  the  output  sensor  array  is  assumed  to  be  of  one 
frequency.  (Recall  the  assumption  made  in  Chapter  3  to  facilitate  double  time 
integration  of  the  acceleration  signals.)  Because  the  output  intensity  is  calculated 
using  this  single  frequency  assumption,  the  input  intensity  is  also  calculated  using  this 
assumption,  to  facilitate  comparison.  Therefore,  the  input  intensity  is  calculated  by 
transforming  both  the  acceleration  and  the  force  signals  into  the  frequency  domain. 
Then  both  signals  are  multiplied  by  Both  signals  are  then  divided  by  (jtoof,  so 
that  the  result  is  to  multiply  both  signals  at  all  frequencies  by  (uVtoo^.  This  has  the 
effect  of  twice  differentiating  the  two  signals  with  a  broadband  assumption,  followed 
by  double  integrating  the  two  signals  with  a  single  frequency  assumption.  The 
acceleration  signal  is  then  divided  by  /to  in  order  to  integrate  to  a  velocity  signal. 
Then,  the  force  and  velocity  signals  are  transformed  to  the  time  domain.  The  two 
time  signals  are  then  multiplied,  and  this  input  intensity  signal  is  transformed  to  the 
frequency  domain,  for  comparison. 

The  output  intensity  is  calculated  on  the  DSP  96002  and  output  via  the  output 
board.  The  output  intensity  is  the  total,  instantaneous  intensity  measured  at  the 
accelerometer  array  and  calculated  as  described  in  Chapter  3.  The  appropriate  DSP 
96002  program  is  run  and  the  intensity  signal  is  input  to  the  analyzer.  The  analyzer 
then  performs  an  FFT  on  the  signal  in  order  to  allow  a  frequency  domain  comparison 
with  the  measured  input  intensity. 


The  two  frequency  representations,  input  and  output  intensity,  are  compared 
for  relative  magnitude  and  phase.  A  discussion  of  the  data  obtained  from  these 
measurements  is  presented  in  Section  6.1. 

5.2.2  Control 

This  section  describes  the  experimental  procedure  involved  in  the  study  of  the 
adaptive  control  of  structural  intensity.  Section  5.2.2. 1  describes  the  data  which  is 
recorded,  and  Section  5.2.2.2  describes  the  procedure  of  data  acquisition. 

5.2.2. 1  Description  of  Control  Data 

The  study  investigates  the  ability  of  the  algorithms  developed  to  control 
intensity.  The  study  also  investigates  the  efrectiveness  of  the  algorithms  relative  to 
the  control  of  acceleration.  The  gauge  of  this  effectiveness  is  an  acceleration  signal 
close  to  the  termination,  referred  to  as  the  downstream  accelerometer.  This  method 
of  evaluation  is  chosen  because,  as  is  noted  in  Section  2.2,  controlling  intensity  should 
greatly  reduce  the  ferfield  vibration.  The  downstream  acceleration  signal  is  obtained 
from  a  PCB  303A13  accelerometer,  and  the  signal  is  conditioned  using  the  PCB 
amplifier  described  above. 

Adaptive  vibration  control  is  performed  at  six  di^rent  frequencies.  Of  these 
six,  three  are  at  resonances  of  the  beam  structure:  140  Hz,  296  Hz,  and  488  Hz.  The 
other  three  frequencies  are  off-resonance:  200  Hz,  340  Hz,  and  520  Hz. 
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Values  of  the  flexural  wavelength  of  vibration  in  the  beam,  K  of  the  nexural 
wavenumber,  k,  of  X/4,  and  of  X/2  are  given  for  each  of  these  six  frequencies  in 
Table  5.1.  If  an  error  sensor  is  at  a  distance  from  the  control  actuator  less  than  V4, 
then  it  is  assumed  to  be  in  the  nearfield  of  the  control  actuator.  At  these  distances, 
the  sensor  will  measure  the  exponentially  decaying  waves  associated  with  the 
nearfield.  If,  however,  the  error  sensor  is  at  a  distance  of  X/2  or  greater  from  the 
control  actuator,  then  it  is  assumed  to  be  in  the  frrfield  of  the  actuator.  The  error 
sensor  is  not  in  the  field  dominated  by  the  exponentially  decaying  waves.  For  the 
beam  in  this  study,  it  is  possible  to  position  the  error  sensor  in  the  nearfield  or  the 
farfield  of  the  control  actuator,  but  positioning  the  error  sensor  in  the  farfield  of  the 
actuator  places  it  in  the  nearfield  of  the  sandbox  termination.  To  avoid  the  confusion 


Table  5.1  Flexural  Wavelength,  Flexural  Wavenumber,  X/4,  and  X/2  of  Beam 
Structure  at  Frequencies  of  Interest. 


Freq.  Hz 

X.  cm 

k  m* 

X/4  cm 

X/2  cm 

140 

116. 

5.43 

29.0 

58.0 

200 

96.9 

6.49 

24.2 

48.5 

296 

79.6 

7.90 

19.9 

39.8 

340 

74.3 

8.46 

18.6 

37.2 

488 

62.0 

10.1 

15.5 

31.0 

520 

60.1 

10.5 

15.0 

30,1 
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that  this  presents,  in  this  thesis  distances  from  the  control  actuator  of  X/2  and  greater 
are  referred  to  as  relatively  far  from  the  control  actuator.  Likewise,  positions  in  the 
nearfield  of  the  control  actuators  are  referred  to  as  relatively  near  to  the  control 
actuators. 

In  an  attempt  to  gain  a  better  understanding  of  the  phenomena  occurring,  the 
position  of  the  control  actuator,  relative  to  the  excitation  shaker  is  varied.  The  two 
positions  of  the  control  actuators  are  referred  to  as  near  to  and  far  from  the  excitation 
shaker.  This  convention  is  used  without  consideration  of  the  nearfield  or  farfleld  of 
the  excitation  shaker. 

There  are  four  hardware  setups  investigated.  These  are  diagrammed  in 
Figure  5.3,  where  X  is  the  distance  from  the  excitation  shaker  to  the  center  error 
accelerometer  and  Y  is  the  distance  from  the  excitation  shaker  to  the  intensity  control 
actuator.  The  distance  from  the  intensity  control  actuator  to  the  center  error 
accelerometer  is,  therefore,  X-Y.  The  flrst  setup  is  confrgured  with  the  control 
shakers  far  from  the  excitation  shaker  (Y=31  cm);  it  also  has  the  error  accelerometer 
array  relatively  far  from  the  control  shakers  (X-Y=60  cm).  The  second  setup  has  the 
control  shakers  positioned  far  from  the  excitation  shaker  (Y=31  cm),  but  the  error 
accelerometer  array  is  relatively  near  to  the  control  shakers  (X-Y=1S  cm).  The  third 
setup  is  configured  with  the  control  shakers  near  to  the  excitation  shaker  (Y=16  cm), 
and  it  has  the  error  accelerometer  array  relatively  far  from  the  control  shakers  (X- 
Ys75  cm).  The  fourth  setup  has  the  control  shakers  positioned  near  to  the  excitation 
shaker  (Y=16  cm),  and  it  has  the  error  accelerometer  array  relatively  near  to  the 


54 


control  shakers  (X-Y=15  cm).  These  set-ups  are  intended  to  investigate  the  influence 
of  farfield  effects  and  nearfield  effects. 


Figure  5.3  Schematic  of  Four  Control  Setups. 

At  each  of  the  twenty-four  setup/firequency  permutations,  seven  control 
strategies  ate  implemented:  (1)  acceleration  control,  (2)  and  (3)  intensity  control 
using  a  force  actuator  only  for  each  of  the  two  gradients  developed  in  Section  4.2.2, 
(4)  and  (5)  intensity  control  using  a  moment  actuator  only  for  each  of  the  two 
gradients,  and  (6)  and  (7)  intensity  control  using  both  force  and  moment  actuators  for 
each  of  the  two  gradients.  This  amounts  to  168  control  scenarios. 
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The  data  recorded  for  each  control  scenario  consists  of  the  control  error  signal 
(speciflc  to  the  control  scenario),  the  structural  intensity  calculated  from  the 
accelerations  measured  from  the  sensor  array,  and  the  downstream  acceleration.  The 
same  data  without  control  is  also  taken  at  each  setup/hrequency  permutation, 
consisting  of  the  acceleration  error  signal,  each  of  the  two  intensity  control  effective 
error  signals,  total  intensity,  and  downstream  acceleration. 

Upon  initial  implementation  of  the  various  intensity  control  schemes,  it  is 
noticed  that  the  error  function  signals  quite  often  contain  high  frequency  noise.  This 
noise  is  seen  at  harmonic  frequencies  of  the  excitation  frequency.  On  certain  tests, 
the  harmonic  noise  is  as  high  or  higher  than  the  signal  at  the  excitation  frequency. 
In  order  to  eliminate  this  high  frequency  noise,  the  error  function  signal  is  passed 
through  a  digital  lowpass  friter.  This  digital  lowpass  filter  is  designed  using  a 
Hamming  window  as  described  in  Oppenheim  and  Schafer  [26:  444-457].  A 

different  filter  is  designed  for  each  frequency  considered. 

5.2.2.2  Description  of  Control  Data  Acquisition 

This  section  describes  the  routine  followed  when  collecting  the  data  described 
in  section  5.2.2. 1.  All  of  the  data  are  recorded  and  analyzed  using  the  HP  3566A 
analyzer  described  above.  For  each  data  signal  recorded,  a  uniform  average  of  50 
power  spectra  is  recorded. 

At  each  setup/frequency  permutation,  the  PCB  conditioning  amplifier  gains  ate 
set  to  ensure  that  the  ±  2.5  V  maximum  range  of  the  input  board  is  not  exceeded. 


56 


The  five  channels  conditioning  the  error  accelerometer  array  are  set  to  the  same  gain. 

After  the  conditioning  amplifier  gains  are  set,  the  data  are  saved  without 
control.  Following  this,  each  of  the  control  scenarios  is  implemented  in  turn.  With 
a  given  scenario  running,  the  intensity  signal,  the  appropriate  error  signal,  and  the 
downstream  acceleration  signal  are  recorded.  A  discussion  of  the  results  obtained 
from  these  measurements  is  given  in  Section  6.2. 


Chapter  6 

DISCUSSION  OF  EXPERIMENTAL  RESULTS 


This  chapter  examines  the  experimental  results  obtained  in  the  study.  Section 
6.1  describes  the  data  obtained  in  the  intensity  measurements.  Section  6.2  describes 
the  adaptive  vibration  control  achieved  through  the  Hltered-x  algorithm. 

6.1  Discussion  of  Intensity  Measurement 

This  section  discusses  the  validity  of  the  intensity  measurements.  The  data 
supporting  these  measurements  are  acquired  in  the  manner  described  in  Section  5.2.1. 

The  relative  magnitudes  of  the  two  intensity  signals  are  discussed  Hrst.  The 
absolute  magnitudes  are  not  discussed  because  the  output  intensity  signal  is  scaled  to 
utilize  the  resolution  of  the  output  board.  The  magnitude  of  the  output  intensity  is 
shown  in  Figure  6.1,  and  the  magnitude  of  the  input  intensity  is  shown  in  Figure  6.2. 
When  one  looks  at  these  two  traces  four  things  require  explanation:  the  noise  levels 
of  the  two  traces,  the  general  positive  slope  of  the  input  intensity  trace,  the  existence 
of  a  dominant  300  Hz  peak  on  the  output  intensity  trace,  and  the  presence  of  a  large 
dc  peak  on  the  input  intensity  trace. 

The  first  issue  to  be  discussed  is  that  the  noise  appears  to  be  higher  on  the 
output  intensity  trace.  This  is  reasonable  because  of  the  greater  number  of  transducers 
involved  in  this  measurement  and  because  of  the  number  of  approximations  and 
calculations  involved  in  the  finite  difference  schemes  used. 
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The  second  issue  for  discussion— the  positive  slope  of  the  input  intensity  trace- 
can  be  attributed  to  the  method  used  to  determine  the  measured  output  intensity.  The 
higher  noise  level  associated  with  the  output  intensity  covers  up  any  slope  in  that 
trace.  Also,  multiplying  the  input  signals  by  (ju>)^  would  introduce  a  slope  of  12 
dB/octave,  which  is  comparable  to  the  slope  appearing  in  the  input  intensity  trace. 

The  third  issue-the  dominant  300  Hz  peak  on  the  output  intensity  ttace-is 
thought  to  be  error  introduced  through  use  of  the  finite  difference  approximations. 

Finally,  the  fourth  issue  for  discussion— the  large  dc  peak  on  the  input  trace-is 
thought  to  be  due  to  the  large  dc  bias  of  the  force  transducer.  With  no  excitation,  the 
force  transducer  gives  a  dc  peak  approximately  75  dB  above  the  noise  level.  This  is 
opposed  to  a  representative  accelerometer  which,  with  no  excitation,  gives  a  dc  peak 
of  approximately  45  dB  above  the  noise  level.  Because  the  output  intensity  is 
measured  only  with  accelerometers  the  dc  bias  would  not  be  as  great. 

The  phases,  calculated  as  arctan(Im/Re),  of  the  two  signals  are  also 
investigated.  The  phase  of  the  output  intensity  is  shown  in  Figure  6.3,  and  the  phase 
of  the  input  intensity  is  shown  in  Figure  6.4.  When  these  phases  are  compared,  one 
can  see  that,  at  the  three  important  frequencies  (0,  300,  and  600  Hz),  the  phases  differ 
by  approximately  180**.  This  is  due  to  the  feet  that  the  positive  direction  for  the  input 
intensity  was  opposite  to  that  used  for  the  output  intensity.  Taking  this  difference  in 
orientation  into  account,  the  phase  differs  by  5*"  or  6°. 

The  output  intensity  is  thought  to  be  a  valid  measure  due  to  this  correlation 
with  input  intensity.  The  feet  that  the  phases  are  almost  matched  and  the  general 


agreement  of  the  magnitude  peaks  support  this  claim.  There  is,  however,  a  great  deal 
of  noise  introduced  into  the  intensity  signal  measured  by  the  accelerometer  array. 

6.2  Discussion  of  Control  Results 

The  results  obtained  from  the  various  control  scheme  measurements  are 
discussed  in  this  section.  There  are  three  items  which  should  be  noted.  The  first  is 
that  for  convenience  in  examining  the  following  frequency  domain  plots,  the  vertical 
divisions  are  placed  at  multiples  of  the  driving  frequency.  The  second  is  that  due  to 
limits  of  the  output  board,  the  intensity  and  error  function  signals  are  scaled  to  avoid 
clipping  while  realizing  the  resolution  potential  of  the  D/A  converter.  This  scaling 
of  the  signals  is  invariant  within  a  given  setup/frequency  permutation.  The  result  is 
that  specifrc  values  are  compared  within  a  setup/frequency  permutation;  however,  only 
relative  reductions  are  to  be  compared  otherwise.  The  third  point  to  note  is  that  for 
convenience  of  discussion,  the  gradients  of  intensity  that  were  developed  in  Section 

4.2.2  are  referred  to  as  the  long  error  function  and  the  short  error  function.  The  error 
function  expressed  in  Eq.  (4.S2)  is  referred  to  as  the  long  error  function: 

e{m)^E!tL*v*F*Et-*Q*kM.  (6.1) 

“o  "o 

The  error  function  expressed  in  Eq.  (4.57)  is  referred  to  as  the  short  error  function: 
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e(/n)=^v+C2.  (6-2) 

There  are  seven  aspects  of  the  control  results  which  are  examined.  Section 
6.2.1  is  a  discussion  of  the  acceleration  control,  which  works  better  with  the  error 
sensor  positioned  relatively  far  from,  rather  than  relatively  near  to,  the  control 
actuator.  Section  6.2.2  discusses  the  general  trends  of  intensity  control  with  the  error 
sensor  positioned  relatively  far  from  the  control  actuator.  Section  6.2.3  examines 
intensity  control  with  the  error  sensor  positioned  relatively  near  to  the  control 
actuator.  Section  6.2.4  is  a  discussion  of  the  general  trend  of  frequency  dependence 
of  the  intensity  control.  Section  6.2.5  is  a  discussion  of  the  issues  associated  with  the 
long  error  function,  while  Section  6.2.6  examines  the  issues  associated  with  the  short 
error  function.  Finally,  Section  6.2.7  discusses  the  trends  in  actuator  configuration 
dependence  in  intensity  control.  It  describes  which  of  the  actuator  confrgurations-- 
force  actuator  only,  moment  actuator  only,  or  force  and  moment  actuators-yields  the 
best  performance. 

The  various  hardware  setups  that  were  examined  in  Section  S.2.2.1  are  given 
in  Table  6.1  for  quick  reference. 

6.2.1  Acceleration  Control 

This  section  examines  the  differences  between  controlling  acceleration  with 
the  error  sensor  positioned  relatively  far  from  and  relatively  near  to  the  control  actua¬ 
tor.  In  Section  2.2  it  is  noted  that  the  control  should  perform  better  with  the  error 


Table  6.1  Hardware  setups  versus  positions  of  error  sensor  and  control  actuator. 


Hardware 
Setup  # 

Error  Sensor  Relative  to 
Control  Actuator 

Control  Actuator  Relative 
to  Excitation  Shaker 

1 

far 

far 

2 

near 

far 

3 

far 

near 

4 

near 

near 

sensor  positioned  relatively  far  from,  rather  than  relatively  near  to,  the  control 
actuator. 

If  one  examines  Figures  6.5  and  6.6,  it  can  be  seen  that  the  control  scheme 
reduces  the  error  acceleration  signal  comparably.  In  each  case  the  error  acceleration 
signal,  at  the  driving  frequency,  was  reduced  by  over  60  dB.  However,  if  one 
examines  Figures  6.7  and  6.8,  this  comparable  reduction  in  the  error  acceleration  does 
not  result  in  comparable  reductions  in  the  downstream  acceleration  signal.  In  Figure 
6.7  it  can  be  seen  that  controlling  acceleration  with  the  error  sensor  positioned 
relatively  far  from  the  control  actuator  results  in  more  than  a  50  dB  drop  in  the 
downstream  acceleration  signal,  while  in  Figure  6.8  it  can  be  seen  that  controlling 
acceleration  with  the  error  sensor  positioned  relatively  near  to  the  control  actuator 
results  in  approximately  a  20  dB  drop  in  the  downstream  acceleration  signal.  These 
results  are  in  agreement  with  the  expectations  discussed  in  Section  2.2. 
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Figure  6.5  Error  Acceleration,  Acceleration  Control,  296  Hz,  Setup  3. 


Figure  6.6  Error  Acceleration,  Acceleration  Control,  296  Hz,  Setup  4. 
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Figure  6.8  Downstream  Acceleration,  Acceleration  Control,  296  Hz,  Setup  4. 
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6.2.2  Intensity  Control  in  the  Farfleld 

This  section  discusses  intensity  control  with  the  error  sensor  positioned 
relatively  far  from  the  control  actuators.  It  is  noted  in  Section  2.2  that  the  intensity 
control  should  perform  as  well  as  the  acceleration  control  with  the  error  sensor 
located  in  the  farfieid. 

If  one  examines  Figure  6.9,  it  can  be  seen  that  the  short  error  function  signal 
is  reduced  by  approximately  50  dB.  One  can  observe  the  roll  off  of  the  error  function 
signal  as  a  result  of  the  lowpass  digital  filter  implemented  as  described  in  Section 
5.2.2.I.  The  intensity  signal,  as  shown  in  Figure  6.10,  is  reduced  as  well.  The  peak 


Figure  6.9  Error  Function,  Short  Error  Intensity  Control,  Force  and  Moment 
Actuators,  488  Hz,  Setup  3. 
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Figure  6.10  Intensity  Signal,  Short  Error  Intensity  Control,  Force  and  Moment 
Actuators,  488  Hz,  Setup  3. 


at  IfOo,  976  Hz,  is  reduced  on  the  order  of  55  dB,  while  the  dc  value  is  reduced 
approximately  25  dB.  The  downstream  acceleration  signal  is  reduced  approximately 
20  dB,  see  Figure  6.11,  through  implementation  of  this  intensity  control  scheme.  This 
is  in  contrast  to  the  more  than  40  dB  reduction  in  the  downstream  acceleration 
achieved  through  the  implementation  of  the  acceleration  control  scheme,  as  shown  in 
Figure  6.12. 

The  fact  that  the  acceleration  control  performs  better  than  the  intensity  control 
when  the  error  sensor  is  positioned  relatively  fiir  from  the  control  actuators  is 
anticipated.  This  is  because  of  the  noise  introduced  into  the  error  function  signal  due 


to  the  numerical  noise  associated  with  the  finite  difference  schemes. 
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Figure  6.11  Downstream  Acceleration,  Short  Error  Intensity  Control,  Force  and 
Moment  Actuators,  488  Hz,  Setup  3. 


Figure  6.12  Downstream  Acceleration,  Acceleration  Control,  488  Hz,  Setup  3. 
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6.23  Intensity  Control  in  the  Nearfieid 

This  section  examines  the  implementation  of  the  intensity  control  with  the 
error  sensor  positioned  relatively  near  to  the  control  actuators.  The  theory  in  Section 
2.2  leads  one  to  expect  that  the  intensity  control  will  outperform  the  acceleration 
control  for  this  conflguration.  The  results  obtained  over  the  range  of  configurations 
are  not  conclusive,  but  seem  to  indicate  that  it  is  possible  to  obtain  improved 
attenuation  using  intensity  control. 

A  representative  example  of  the  adaptive  control  of  the  short  error  function 
signal  is  given  in  Figure  6.13.  The  error  function  is  controlled  well  in  most  cases. 
The  reduction  is  usually  on  the  order  of  40  to  50  dB.  In  this  case,  the  corresponding 
intensity  signal  experiences  a  significant  reduction.  As  shown  in  Figure  6.14,  the 
intensity  signal  is  reduced  approximately  30  dB  at  the  frequency  20)^,  and  it  is 
reduced  approximately  30  dB  at  the  dc  value.  The  resulting  reduction  of  the 
downstream  acceleration  signal,  shown  in  Figure  6.15,  is  on  the  order  of  17  dB.  This 
is  compared  with  the  downstream  acceleration  reduction  achieved  with  acceleration 
control,  shown  in  Figure  6.16,  of  approximately  10  dB.  Another  example  of  intensity 
control  performing  better  than  acceleration  control  is  shown  next.  In  Figure  6.17,  the 
downstream  acceleration  signal  resulting  from  intensity  control  at  340  Hz  is  shown. 
The  reduction  in  downstream  acceleration  is  approximately  15  dB.  This  is  compared 
to  the  reduction  in  downstream  acceleration  resulting  from  acceleration  control,  which 
is  shown  in  Figure  6.18.  The  reduction  is  approximately  5  dB. 


Figure  6.14  Intensity  Signal,  Short  Error  Intensity  Control,  Moment  Actuator,  488 
Hz,  Setup  2. 
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Figure  6.15  Downstream  Acceleration,  Short  Error  Intensity  Control,  Moment 
Actuator,  488  Hz,  Setup  2. 


Figure  6.16  Downstream  Acceleration,  Acceleration  Control,  488  Hz,  Setup  2. 
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Figure  6.17  Downstream  Acceleration,  Short  Error  Intensity  Control,  Force  and 
Moment  Actuators,  340  Hz,  Setup  2. 


Figure  6.18  Downstream  Acceleration,  Acceleration  Control,  340  Hz,  Setup  2. 
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It  should  be  noted  that  the  cases  where  the  intensity  control  outperforms  the 
acceleration  control  are  few.  In  the  majority  of  the  tests  run,  the  intensity  control 
falls  short  of  the  acceleration  control.  In  performing  the  tests,  several  trends  are 
identified  which  may  contribute  to  these  difficulties.  Some  of  these  trends  and 
problems  are  discussed  in  the  following  sections. 

6.2.4  Frequency  Dependence  of  Intensity  Control 

The  first  problem  associated  with  intensity  control  that  is  discussed  is  that  of 
frequency  dependence.  The  short  error  intensity  control  scheme  with  a  force  actuator 
only  is  examined  at  all  frequencies.  The  control  actuator  is  positioned  near  to  the 
excitation  shaker  and  the  error  sensor  is  positioned  relatively  near  to  the  control 
actuators. 

The  control  scheme  works  better  at  the  higher  frequencies.  This  seems  to  be 
true  because  the  error  function  has  a  higher  signal-to-noise  ratio  at  these  frequencies. 
In  Figure  6.19,  it  can  be  seen  that  although  the  control  reduced  the  error  function 
signal  to  the  noise  floor,  at  the  140  Hz  excitation  frequency,  this  amounts  to  only  a 
25  dB  reduction  in  the  error  function  signal.  This  results  in  a  small  reduction  in  the 
downstream  acceleration  signal,  as  shown  in  Figure  6.20.  Likewise,  as  shown  in 
Figure  6.21  which  displays  control  with  an  excitation  frequency  of  340  Hz,  a  signal- 
to-noise  ratio  for  the  error  function  signal  of  about  35  dB  re  Its  in  an  increase  in  the 
downstream  acceleration  signal,  shown  in  Figure  6.22.  However,  when  the  signal-to- 
noise  ratio  of  the  error  function  is  on  the  order  of  45  to  50  dB,  as  is  the  case  for 
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Figure  6.19  Error  Function,  Short  Error  Intensity  Control,  Force  Actuator,  140  Hz, 
Setup  4. 


Figure  6.20  Downstream  Acceleration,  Short  Error  Intensity  Control,  Force  Achiator, 
140  Hz,  Setup  4. 
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Figure  6.21  Error  Function,  Short  Error  Intensity  Control,  Force  Actuator,  340  Hz, 
Setup  4. 


Figure  6.22  Downstream  Acceleration,  Short  Error  Intensity  Control,  Force  Actuator, 
340  Hz,  Setup  4. 
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excitation  frequencies  296,  488,  and  520  Hz  (see  Figures  6.23,  6.24,  and  6.25, 
respectively),  the  resulting  reduction  in  the  downstream  acceleration  signal  is  on  the 
order  of  10-20  dB,  as  shown  in  Figures  6.26,  6.27,  and  6.28.  It  should  be  noted  that 
at  200  Hz,  the  control  scheme  failed  to  converge.  This  strengthens  the  claim  that  the 
scheme  works  better  at  higher  frequencies. 

A  possible  explanation  for  this  frequency  dependence  concerns  the  accelerome¬ 
ter  spacing  in  terms  of  the  flexural  wavelength,  k.  The  wavelengths  for  the  various 
frequencies  of  interest  are  given  in  Table  5.1.  The  spacing  between  accelerometers 
is  4  cm.  At  the  highest  frequency,  520  Hz,  this  corresponds  to  about  0.067  k.  At  the 
lowest  frequency,  140  Hz,  this  accelerometer  spacing  corresponds  to  about  0.034  k. 
The  reason  the  accelerometers  should  be  spaced  farther  apart,  relative  to  a  wave¬ 
length,  is  that  in  the  spatial  flnite  differencing  schemes  neighboring  acceleration 
values  are  subtracted,  in  order  to  approximate  the  spatial  derivatives.  If  these  signals 
are  almost  equal,  more  error  will  be  introduced,  with  a  finite  bit  accuracy  processor, 
than  if  the  signals  are  not  almost  equal.  Positioning  the  accelerometers  close,  relative 
to  a  wavelength,  ensures  the  neighboring  values  will  be  almost  equal. 

Another  possible  explanation  for  the  frequency  dependence  concerns  the 
sampling  frequency  versus  the  driving  frequency.  This  discussion  is  analogous  to  the 
one  for  geometrical  spacing  given  previously.  The  successive,  in  time,  signals  are 
subtracted  to  give  an  approximation  of  the  time  derivatives.  The  closer  these  signals 
are  to  being  equal,  the  more  processing  error  will  be  introduced  into  the  derivative 
approximation.  Sampling  the  signals  closer,  relative  to  driving  frequency. 
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Figure  6.23  Error  Function,  Short  Error  Intensity  Control,  Force  Actuator,  296  Hz, 
Setup  4. 


Figure  6.24  Error  Function,  Short  Error  Intensity  Control,  Force  Actuator,  488  Hz, 
Setup  4. 
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Figure  6.25  Error  Function,  Short  Error  Intensity  Control,  Force  Actuator,  520  Hz, 
Setup  4. 


Figure  6.26  Downstream  Acceleration,  Short  Error  Intensity  Control,  Force  Actuator, 
296  Hz,  Setup  4. 
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Figure  6.27  Downstream  Acceleration,  Short  Error  Intensity  Control,  Force  Actuator, 
488  Hz,  Setup  4. 


Figure  6.28  Downstream  Acceleration,  Short  Error  Intensity  Control,  Force  Actuator, 
520  Hz,  Setup  4. 


ensures  that  the  successive  samples  are  closer  to  being  equal. 

A  third  explanation  for  this  frequency  dependence  relates  the  beam  length  to 
the  wavelength.  One  should  recall  that  a  +x  travelling  wave  is  assumed  in  the 
development  of  the  intensity  error  functions.  The  wavelength  at  lower  frequencies 
approaches  the  length  of  the  beam.  Thus,  the  error  sensor  is  in  a  different  wave  field 
at  the  lower  frequencies  than  at  the  higher  frequencies.  The  more  the  Held  is 
dominated  by  waves  other  that  4-jc  travelling  waves,  the  more  the  error  function  is  no 
longer  valid,  and  the  less  the  calculated  error  function  relates  to  the  true  gradient  of 
intensity. 

6.2.5  Long  Error  Function 

The  next  problem  associated  with  intensity  control  concerns  the  long  error 
function.  A  reduction  in  the  long  error  function  signal  used  does  not  necessarily 
result  in  a  reduction  in  the  intensity  signal. 

As  is  shown  in  Figure  6.29,  the  given  control  scheme  succeeded  in  reducing 
the  error  function  signal  by  approximately  55  dB.  However,  the  intensity  signal, 
shown  in  Figure  6.30,  increased  to  the  point  that  it  exceeded  the  ±8.19  V  limit  of  the 
output  board.  It  should  be  noted  that  because  the  intensity  signal  is  clipped,  the 
frequency  plot  of  Figure  6.30  is  not  a  true  representation  of  the  signal,  but  it  serves 
to  express  the  fact  that  the  intensity  signal  increased. 


Figure  6.30  Intensity  Signal,  Long  Error  Intensity  Control,  Force  Actuator,  520  Hz, 
Setup  2. 
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This  apparent  lack  of  correlation  between  the  long  error  function  and  the 
intensity  is  not  universal.  As  is  shown  in  Figures  6.31  and  6.32,  controlling  the  long 
error  signal  sometimes  reduces  the  intensity. 

One  reason  for  this  correlation/lack  of  correlation  concerns  the  +x  travelling 
wave  assumption  made  in  estimating  the  gradient.  If  the  wave  field  is  dominated  by 
another  type  of  wave,  the  expression  for  the  gradient  does  not  necessarily  correlate 
with  the  true  gradient,  and,  therefore,  it  does  not  correlate  with  the  intensity. 


Figure  6.31  Error  Function,  Lx}ng  Error  Intensity  Control,  Force  Actuator,  520  Hz, 
Setup  3. 
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Figure  6.32  Intensity  Signal,  Long  Error  Intensity  Control,  Force  Actuator,  520  Hz, 
Setup  3. 


6.X6  Short  Error  Function 

It  is  found  that,  although  the  short  error  function  seems  to  correlate  with  the 
intensity,  there  are  still  problems  with  this  error  function.  One  problem  is  that  the 
error  signal  is  too  noisy  at  low  frequencies.  As  can  be  seen  in  Figure  6.33,  the 
signal-to-noise  ratio  of  the  error  signal  at  200  Hz  for  setup  1  is  less  than  30  dB.  This 
noise  is  thought  to  be  numerical  noise  resulting  from  the  finite  difference  schemes. 
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Figure  6.33  Short  Error  Function,  No  Control,  200  Hz,  Setup  1. 


6.2.7  Actuator  Configuration  Dependence 

The  Enal  issue  concerning  intensity  control  to  be  examined  is  that  of  control 
actuator  configuration  dependence.  There  are  no  strong  trends  of  one  actuator 
configuration  outperforming  another  to  any  great  degree,  but  the  trend  appears  to  be 
that  the  moment  actuator  conHguration  works  slightly  better  than  the  force  and 
moment  actuators  configuration,  which  works  slightly  better  than  the  force  actuator 
configuration. 

At  the  frequency  examined,  the  short  error  function  signal  for  each  of  the  three 
configurations  is  reduced  by  similar  amounts.  A  representative  case  of  this  is  shown 
in  Figure  6.34;  the  error  function  in  all  cases  is  reduced  on  the  order  of  45  dB. 


However,  the  difference  is  found  in  the  traces  of  downstream  acceleration.  The  force 
actuator  configuration,  shown  in  Figure  6.35,  reduces  the  downstream  acceleration 
about  13  dB.  The  moment  actuator  configuration,  shown  in  Figure  6.36,  reduces  the 
downstream  acceleration  by  approximately  17  dB.  The  force  and  moment  actuators 
configuration,  shown  in  Figure  6.37,  reduces  the  downstream  acceleration  slightly 
less—about  15  dB. 
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Figure  6.34  Error  Function,  Short  Error  Intensity  Control,  Force  Actuator,  340  Hz, 
Setup  2. 
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Figure  6.37  Downstream  Acceleration,  Short  Error  Intensity  Control,  Force  and 
Moment  Actuators,  340  Hz,  Setup  2. 


Chapter  7 

CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  chapter.  Section  7.1  examines  the  conclusions  of  the  research,  and 
Section  7.2  recommends  several  possible  areas  for  future  research. 

7.1  Conclusions 

The  goals  of  this  thesis  were  to  develop,  implement,  and  evaluate  several 
algorithms  for  the  adaptive  vibration  control  of  the  total,  instantaneous  structural 
intensity  in  a  beam.  A  measurement  algorithm  based  on  Hnite  difference  schemes 
was  developed  for  calculating  the  necessary  partial  derivatives  in  both  space  and  time. 
A  control  algorithm  based  on  the  filtered-x  algorithm  was  developed.  The  intensity 
control  algorithm  used  one  of  two  derived  error  functions  and  one  of  three  actuator 
configurations  to  implement  the  control  of  structural  intensity.  The  various  intensity 
control  schemes  were  evaluated  relative  to  acceleration  control  by  means  of  a 
downstream  acceleration  signal.  Evaluation  of  the  results  from  these  tests  led  to 
several  conclusions. 

The  first  conclusion  is  that  acceleration  control  was  more  effective  with  the 
error  sensor  in  the  faifield  than  in  the  nearfield,  which  was  predicted  in  Section  2.2. 

Furthermore,  with  the  error  sensors  positioned  relatively  far  from  the  control 
actuators,  controlling  acceleration  was  more  effective  than  controlling  intensity,  using 
the  algorithms  developed  in  this  study.  Although  the  theory  predicts  that  the  two 
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schemes  will  perform  equally  well,  noise  in  the  error  signal  degrades  the  effectiveness 
of  the  control  algorithm,  when  controlling  intensity. 

In  addition,  with  the  error  sensors  positioned  relatively  near  to  the  control 
actuators,  the  attenuation  achieved  by  controlling  intensity  was  comparable  to  or 
greater  than  that  achieved  by  control  ting  acceleration.  It  was  predicted  that  this 
attenuation  should  always  be  greater  when  controlling  intensity  than  when  controlling 
acceleration.  The  reasons  for  this  discrepancy  include  frequency  dependence  of  the 
control  schemes  and  numerical  noise  in  the  computation  of  the  error  function  signals. 

The  performance  of  the  intensity  control  schemes  is  frequency  dependent  due 
to  the  frequency  dependence  of  the  numerical  noise  introduced  by  the  finite  difference 
schemes.  In  addition,  the  beam  length  relative  to  the  flexural  wavelength  of  vibration 
at  low  frequencies  places  the  error  sensor  in  a  wave  field  dominated  by  waves  other 
than  the  assumed  -t-x  travelling  waves. 

The  performance  of  the  intensity  control  schemes  is  strongly  dependent  upon 
the  choice  of  error  function.  The  long  error  function  frequently  does  not  give  a  true 
representation  of  the  gradient  of  intensity.  This  is  due  to  the  assumption  that  the 
wave  field  is  dominated  by  +x  travelling  waves  which  is  not  valid  in  the  nearfield. 

Signal  noise  was  a  problem  at  all  frequencies,  with  both  gradient  approxima¬ 
tions.  This  noise  consists  of  electronic  noise  from  the  accelerometers  and  numerical 
noise  from  the  finite  difference  techniques. 

The  performance  of  the  intensity  control  schemes  does  not  seem  to  depend 
upon  the  actuator  configuration.  The  three  configurations,  force  only,  moment  only. 
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and  force  and  moment,  all  result  in  similar  reductions  in  both  the  error  function  signal 
and  the  downstream  acceleration  signal. 

7.2  Recommendations  for  Further  Research 

Based  upon  the  results  of  this  study  and  the  above  conclusions,  several 
recommendations  for  future  research  are  made.  These  should  include  a  more  in-depth 
study  of  the  intensity  measurement,  including  the  gradient  of  intensity,  developing  the 
control  algorithms  to  be  able  to  discriminate  between  active  and  reactive  intensity,  and 
increasing  the  complexity  of  the  structure  being  controlled. 

A  different  method  for  calculating  the  intensity  will  be  needed  if  control  of 
broadband  vibrations  is  to  be  achieved.  The  method  described  in  this  thesis  used  a 
multiplication  by  to  achieve  double  integration  in  time,  which  assumes  single 
frequency  excitation  at  a  known  ftequency,  tO),.  One  alternative  would  be  to 
investigate  using  piezoelectric  strain  gauges.  These  strain  gauges  would  allow  a 
measure  of  strain,  which  may  lead  to  an  improved  intensity  estimate  capable  of 
making  broadband  intensity  measurements.  Another  alternative  would  be  to  pass  the 
acceleration  signals  through  linear  analog  integrators,  designed  to  be  effective  over 
the  frequency  range  of  interest.  Using  transducer  combinations,  such  as  both  strain 
gauges  and  accelerometers,  should  also  be  considered. 

It  should  be  noted  that  any  multiple  frequency  study  would  necessitate  a 
broadband  calibration  of  the  transducers.  Currently,  the  control  schemes  use  a  single 
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calibration  constant  for  each  accelerometer  at  the  excitation  frequency.  Perhaps  a 
frequency  domain  calibration  could  be  performed  using  a  FIR  filter  approximation. 

Another  quantity  that  warrants  further  investigation  is  the  estimate  of  the 
gradient  of  the  intensity.  At  times,  a  substantial  reduction  in  the  long  error  signal 
used  in  this  study  did  not  result  in  any  reduction  in  the  intensity.  A  reevaluation  of 
the  assumptions  will  be  necessary  in  order  to  hopefully  develop  a  more  valid 
expression  for  the  gradient  of  intensity. 

An  additional  area  to  investigate  is  the  separation  of  the  active  and  the  reactive 
intensity.  Because  the  active  intensity  represents  the  energy  which  propagates  to  the 
farfield,  it  would  often  be  desirable  to  control  just  the  active  intensity. 

A  final  recommendation  is  to  extend  the  method  of  controlling  intensity  to 
more  complex  structures.  For  one-dimensional  structures,  this  could  include  different 
terminations,  terminations  at  both  ends,  and  coupled  beams  with  various  joints.  The 
study  could  also  be  implemented  on  two-dimensional  structures,  such  as  plates  and 


shells. 
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Table  A.1  Calibration  Constants  for  Error  Array. 


Freq. 

(Hz) 

Ser.# 

1136 

Ser.# 

1140 

Ser.# 

1141 

Ser.# 

1142 

Ser.# 

1156 

120 

1.0000 

0.8977 

0.9270 

1.850 

1.229 

140 

1.0000 

0.8747 

0.9072 

1.662 

1.245 

160 

1.0000 

0.8243 

0.8860 

1.493 

1.261 

180 

1.0000 

0.8276 

0.8638 

1.371 

1.274 

200 

1.0000 

0.8050 

0.8412 

1.285 

1.281 

220 

1.0000 

0.7842 

0.8187 

1.217 

1.279 

240 

1.0000 

0.7656 

0.7969 

1.156 

1.266 

1 

1.0000 

0.7494 

0.7758 

1.094 

1.241 

1  280 

1.0000 

0.7363 

0.7560 

1.037 

1.209 

300 

1.0000 

0.7260 

0.7379 

0.9897 

1.175 

320 

1.0000 

0.7186 

0.7218 

0.9545 

1.138 

1 

1.0000 

0.7143 

0.7080 

0.9295 

1.107  1 

1  ^ 

1.0000 

0.7120 

0.6959 

0.9107 

1.080  1 

1  380 

1.0000 

0.7113 

0.6858 

0.8925 

1.060  1 

1  ^ 

1.0000 

0.7124 

0.6776 

0.8724 

1.044  1 

420 

1.0000 

0.7141 

0.6705 

0.8502 

1.034  1 

440 

1.0000 

0.7164 

0.6645 

0.8275 

1.028 

460 

1.0000 

0.7189 

0.6591 

0.8044 

1.026 

480 

1.0000 

0.7217 

0.6550 

0.7792 

1.027 

500 

1.0000 

0.7238 

0.6510 

0.7433 

1.031 

520 

1.0000 

0.7261 

0.6474 

0.6930 

1.038 

